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Requirements  for  the  Degree  of  Doctor  of  Philosophy 

SCHEDULING  PROBLEMS  IN  PROJECT  MANAGEMENT 

By 

Oya  Icmeli 
December  1992 
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Major  Department:  Decision  and  Information  Sciences 

In  this  dissertation  we  study  two  new  scheduling  problems 
in  project  management.  The  first  one  is  called  the  Resource 
Constrained  Project  Scheduling  Problem  with  Discounted  Cash 
flows  (RCPSPDC) . Exact  and  heuristic  solution  procedures  for 
this  problem  are  developed.  A Tabu  Search  procedure  is 
proposed  (Tabu-S)  as  a heuristic  solution  technique  for  this 
problem.  The  procedure  is,  then,  modified  to  invoke  a long 
term  memory  function  (Tabu-L) . Both  procedures  are  tested  on 
50  problems  derived  from  Patterson's  data  set.  Solutions 
produced  by  these  procedures  are  compared  to  upper  bounds 
obtained  from  a linear  programming  relaxation  of  RCPSPDC  which 
is  strengthened  by  valid  cuts.  In  general.  Tabu  Search 
successfully  produces  near-optimal  solutions  with  reasonable 
computational  effort. 
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A Branch  and  Bound  solution  procedure  is  developed  to 
obtain  an  optimal  solution  for  RCPSPDC.  The  upper  bounds  at 
the  nodes  of  the  branch  and  bound  tree  are  obtained  by  solving 
payment  scheduling  problems.  It  is  shown  that  this  procedure 
finds  an  optimal  solution  to  RCPSPDC  in  a finite  number  of 
steps.  The  computational  results  on  50  test  problems 
demonstrates  that  the  solutions  are  obtained  within  0.5% 
optimality  with  reasonable  computational  effort. 

The  second  scheduling  problem  that  we  examine  is  the 
Resource  Constrained  Time/Cost  Tradeoff  Project  Scheduling 
Problem  with  Discounted  Cash  flows  (RCTCTPSPDC) . To  our 
knowledge,  this  problem  is  first  introduced  to  project 
scheduling  literature  in  this  dissertation.  For  this  problem 
a bounding  algorithm  and  a heuristic  solution  procedure  are 
developed.  The  heuristic  procedure  performs  a single  forward 
pass  through  the  network  starting  with  the  first  activity  and 
schedules  the  activities  according  to  their  net  marginal  gain 
values  until  the  last  activity  is  scheduled.  The  procedure  is 
tested  on  50  test  problems  and  is  very  efficient  in  that  each 
problem  required  less  than  0.5  second  CPU  time  on  the  IBM  3090 
computer.  Furthermore,  a parametric  analysis  was  done  to  show 
the  effect  of  problem  structure  on  the  results. 

Finally,  a subgradient  algorithm  is  advanced  to  solve  the 
RCTCTPSPDC.  This  algorithm  employs  the  Generalized  Benders 
Decomposition  technique  to  solve  the  subproblems  obtained  when 
the  lagrangian  multipliers  are  fixed.  The  near  optimal 
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CHAPTER  1 
INTRODUCTION 


Project  management  is  known  as  a specialized  branch  of 
management  which  has  evolved  to  plan,  schedule  and  control 
complex  activities  to  achieve  performance,  cost  and  time 
objectives  for  a given  scope  of  work.  Project  planning  starts 
with  preliminary  coordination  where  the  various  parties 
involved  in  the  project  get  together  and  determine  the  project 
objectives.  To  achieve  those  objectives,  the  project  budget 
and  the  detailed  description  of  the  activities  that  must  be 
undertaken  are  determined  [20].  Plans  are  also  developed  to 
deal  with  the  project  termination  and  postcontrol  of  the 
project.  The  planning  stage  is  followed  by  the  implementation 
of  these  plans  which  is  recognized  as  scheduling.  At  this 
stage,  the  starting  and  finishing  times  of  activities,  the 
critical  path  and  float  times  associated  with  the  project  are 
identified.  Furthermore,  adjustments  on  the  schedule  are  made 
in  order  to  allocate  or  smooth  the  use  of  resources,  if  the 
resources  are  scarce.  A number  of  techniques  have  been 
developed  to  assist  in  scheduling  of  projects.  These 
techniques  usually  use  network  models  as  an  aid  in  the 
scheduling.  The  two  popular  methods  are  the  Critical  Path 
Method  (CPM)  and  the  Program  Evaluation  and  Review  Technique 
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(PERT) . If  activity  durations  are  known  with  certainty,  CPM 
can  be  used  to  determine  the  length  of  time  required  to 
complete  a project.  The  technique  CPM  can  also  be  used  to 
determine  how  long  each  activity  can  be  delayed  without 
delaying  the  project  completion  time.  The  activities  which 
can  not  be  delayed  without  delaying  the  project  completion 
time  are  recognized  as  the  critical  activities  and  require 
more  attention  from  the  project  manager.  The  technique  was 
developed  in  the  late  1950s  by  researchers  at  Dupont  and 
Sperry  Rand  [39,43]. 

If  the  activity  durations  are  not  known  with  certainty, 
PERT  can  be  used  to  estimate  the  expected  activity  durations 
and  the  probability  that  the  project  will  be  completed  by  a 
given  deadline.  PERT  was  also  developed  in  the  late  1950s  by 
consultants  working  on  the  development  of  the  Polaris  Missile. 
PERT  and  CPM  are  given  a major  share  of  the  credit  for  the 
growing  interest  in  project  management. 

The  last  stage  of  project  management  is  to  control  the 
project.  Information  is  collected  about  project  performance, 
cost  and  time  and  compared  with  the  desired  level  and  action 
taken  if  actual  achievement  differs  from  the  desired 
significantly  [51,52]. 

In  this  dissertation,  we  focus  on  scheduling  problems  in 
project  management.  In  the  following  section,  the  three 
fundamental  project  scheduling  problems,  the  resource 
constrained  project  scheduling  problem  (RCPSP) , the  time/cost 
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tradeoff  problem  (TCTP) , and  the  payment  scheduling  problem 
(PSP)  are  discussed.  Moreover,  new  problems  which  are 
combinations  of  two  or  more  of  these  three  classes  are 
introduced. 

In  the  late  1950s  the  development  of  PERT  and  CPM 
techniques  allowed  projects  to  be  portrayed  by  network 
diagrams  where  jobs  or  activities  are  represented  by  arcs, 
events  are  represented  by  nodes,  and  the  interrelations 
between  the  jobs  or  activities  are  defined  by  the  network 
structure.  However,  project  scheduling  with  these  techniques 
deals  only  with  the  time  aspect  without  consideration  of 
resource  restrictions  and/or  cash  flows.  Yet,  in  many  real 
life  situations,  delays  in  the  execution  time  of  certain 
activities  occur  when  resources  required  by  these  activities 
are  not  available  in  sufficient  quantities  during  the  time 
interval  when  these  activities  are  scheduled  to  take  place. 
This  particular  problem  is  known  as  "the  resource-constrained 
project  scheduling  problem"  in  the  literature. 

The  project  scheduling  problem  referred  to  as  "the 
time/cost  tradeoff  problem"  in  the  literature  involves  the 
reduction  of  certain  activity  durations  by  allocating  more 
resources  to  these  activities.  This  generally  results  in 
higher  direct  activity  costs  and  shorter  project  durations. 

Another  practical  project  scheduling  problem  considers 
the  situation  where  there  are  cash  inflows  and  outflows 
throughout  the  life  of  a project  and  the  project  manager  tries 
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to  schedule  the  project  activities  in  such  a way  that  the  net 
present  value  of  the  cash  flows  is  maximized.  This  particular 
problem  is  known  as  "the  payment  scheduling  problem"  in  the 
literature. 

Research  on  these  three  scheduling  problems  has  led  to 
the  emergence  of  new  scheduling  techniques.  These  techniques 
often  employ  CPM  and/or  PERT  as  a subroutine  in  trying  to 
attain  a solution  for  the  problem  under  consideration. 

The  combination  of  RCPSP  and  PSP  has  led  to  the 
definition  of  a project  scheduling  problem  recognized  as 
Resource  Constrained  Project  Scheduling  Problem  with 
Discounted  Cash  Flows  (RCPSPDC) . In  many  real  life  projects, 
project  completion  requires  the  use  of  various  resources  whose 
limited  availability  may  extend  the  project  completion  time. 
Therefore,  the  project  manager  is  not  only  concerned  about 
allocation  of  resources  but  also  the  time  value  of  money.  In 
response  to  this  concern  the  RCPSPDC  was  formulated.  The 
problem  RCPSPDC  involves  scheduling  the  project  activities 
with  cash  inflows  and  outflows  in  such  a way  that  the 
precedence  and  resource  restrictions  are  satisfied  with  the 
objective  of  maximizing  the  net  present  value  of  the  cash 
flows. 

The  Resource  Constrained  Time/Cost  Tradeoff  Project 
Scheduling  Problem  (RCTCTP)  which  is  the  combination  of  RCPSP 
and  TCTP , concerns  the  scheduling  of  project  activities  with 
resource  and  precedence  constraints  in  such  a way  that  the 
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project  completion  time  is  minimized  where  the  duration  of 
each  activity  may  vary  according  to  the  level  of  resources 
allocated.  Finally,  the  combination  of  PSP  and  TCTP  is  known 
as  the  Time/Cost  Tradeoff  Project  Scheduling  Problem  with 
Discounted  Cash  Flows  (TCTPDF) . This  problem  involves 
determining  the  activity  durations  and  activity  starting  times 
in  such  a way  that  the  precedence  relations  are  met  and  the 
net  present  value  of  cash  flows  are  maximized. 

The  three  fundamental  scheduling  classes,  RCPSP,  TCTP, 
PSP,  can  also  be  integrated  to  define  the  project  scheduling 
problem  which  can  be  named  as  the  Resource  Constrained 
Time/Cost  Tradeoff  problem  with  Discounted  Cash  Flows 
(RCTCTPDC)  . To  our  knowledge,  there  is  no  work  in  the 
literature  on  this  problem.  This  problem  involves  determining 
the  duration  and  the  starting  time  of  each  project  activity  in 
such  a way  that  the  precedence  and  resource  restrictions  are 
met  and  the  net  present  value  of  the  cash  flows  is  maximized. 

The  following  figure  (Figure  1)  summarizes  the  relation 
among  the  project  scheduling  problems. 
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Figure  1.  Relations  among  the  project  scheduling 
problems. 


CHAPTER  2 
LITERATURE  REVIEW 


In  this  chapter,  a unifying  survey  of  the  literature  on 
the  three  fundamental  project  scheduling  problems  is 
presented:  the  resource  constrained  project  scheduling 
problem  (RCPSP) , the  time/cost  tradeoff  problem  (TCTP) , and 
the  payment  scheduling  problem  (PSP) . The  problems  TCTP  and 
PSP  are  linear  programs  which  can  be  solved  with  relative 
ease.  The  problem  RCPSP  is  0-1  Integer  linear  program  for 
which  efficient  special-purpose  branch  and  bound  technigues 
exist.  The  first  complete  survey  in  this  area  was  done  by 
Davis  [11]  in  1973.  The  primary  focus  of  his  survey  was  the 
resource  constrained  project  scheduling  problem.  This  survey 
focuses  on  the  research  since  1973.  In  addition,  the  project 
scheduling  problems  that  result  from  combining  two  or  more  of 
the  three  fundamental  problems  are  discussed.  Inspite  of 
their  practical  relevance,  very  little  work  has  been  done  on 
these  problems  to  date.  The  future  of  the  project  scheduling 
literature  appears  to  be  developing  in  the  direction  of 
combining  the  fundamental  problems  and  developing  efficient 
exact  and  heuristic  methods  for  the  resulting  problems. 
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2.1  Resource  Constrained  Project  Scheduling  Problem 
Since  the  1960s  most  work  in  project  scheduling  has 
focused  on  developing  solution  techniques  for  RCPSP.  Numerous 
heuristics  and  many  exact  solution  techniques  have  emerged. 
The  first  complete  survey  of  this  area  was  performed  by  Davis 
[11].  In  his  paper,  Davis  categorized  the  resource  allocation 
problems  into  three  types:  time/cost  tradeoff  problems, 

problems  in  which  resource  demands  are  leveled,  and  project 
scheduling  problems  with  fixed  resource  limits.  The  primary 
concern  of  the  survey  was  on  the  last  type  of  problems,  those 
which  involve  scheduling  under  fixed  resource  limits.  In 
addition,  Davis  remarked  on  the  strong  similarities  that  exist 
between  the  project  scheduling  problems  and  job  shop 
sequencing  and  assembly  line  balancing  problems.  The 
correspondence  between  project  scheduling  and  assembly  line 
balancing  problems  was  first  detailed  by  Moodie  and  Mandeville 
[53]  and  summarized  by  Davis  [11]  as  follows: 


Assembly  Line  Balancing 

Work  elements 
Work  element  times 

Work  stations 
Cycle  Time 


Project  Network 
Scheduling 
Activities 
Activity  resource 
requirements 
Days 

Maximum  available 
units  of  resource 


Figure  2.  The  correspondence  between  the  project  scheduling 
and  assembly  line  balancing  problems 
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The  conceptual  similarities  between  the  project 
scheduling  problem  and  the  job  shop  problem  were  given  by  an 
example.  Davis'  paper  is  an  excellent  survey  of  the  pre-1973 
project  scheduling  literature.  Consequently,  this  survey  will 
concentrate  on  the  research  findings  since  1973. 

As  Davis  remarked,  although  the  formulation  of  RCPSP  as  a 
0-1  integer  linear  programming  problem  is  common,  solving  it 
directly  is  impractical  and  unattractive  due  to  its  excessive 
computational  requirements.  In  1974,  Patterson  and  Huber  [58] 
formulated  the  resource-constrained  single  project  scheduling 
problem  as  a 0-1  integer  linear  program.  Their  formulation 
adapted  a model  for  multiple  project  scheduling  previously 
developed  by  Pritsker,  Walters,  and  Wolfe  [62].  Rather  than 
solving  their  formulation  directly,  Patterson  and  Huber 
suggested  two  exact  bounding  algorithms,  a minimum  bounding 
algorithm  and  a maximum  bounding  algorithm,  that  examine  the 
feasibility  of  a series  of  0-1  programming  problems.  The 
minimum  bounding  algorithm  starts  with  an  infeasible  yet 
superoptimal  project  completion  time  and  successively 
increases  the  completion  time  by  one  unit  until  an  optimal 
(i.e.,  feasible)  solution  is  obtained.  Conversely,  the 
maximum  bounding  algorithm  starts  with  a feasible  yet 
suboptimal  project  completion  time  and  iteratively  decreases 
the  completion  time  until  an  infeasible  solution  is  found.  In 
addition,  the  starting  solutions  of  the  two  bounding 
algorithms  were  used  to  define  an  interval  for  the  project 
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completion  time  upon  which  Patterson  and  Huber  performed  a 
binary  search.  Finally,  they  compared  their  approach  with  a 
general-purpose  commercial  Integer  Programming  code  on  a set 
of  problems  from  the  literature,  and  demonstrated  that  their 
scheme  required  less  computational  effort  than  the  direct 
solution  procedure. 

Later,  in  1978,  Stinson,  Davis  and  Khumawala  [73] 
formulated  the  multiple  resource-constrained  scheduling 
problem  as  an  integer  programming  problem  and  advanced  a 
branch  and  bound  algorithm  for  solving  it.  The  algorithm  they 
developed  is  similar  to  Johnson's  [38]  branch  and  bound 
algorithm  with  differences  in  the  node  selection  heuristics 
employed  and  the  number  of  resources  handled  (Johnson's 
algorithm  allows  for  a single  resource) . In  their  algorithm, 
branching  corresponds  to  creating  new  partial  feasible 
schedules  from  given  partial  feasible  schedules.  In  an  effort 
to  reduce  the  size  of  the  branch  and  bound  tree,  they 
suggested  two  pruning  operations,  dominance  and  lower 
bounding.  Node  selection  heuristics  were  employed  along  with 
pruning  operations  to  control  the  growth  of  the  tree.  Instead 
of  using  a single  node-selection  rule,  they  employed  a 
decision  vector  consisting  of  several  partial  schedule 
attributes  (an  idea  similar  to  that  of  Ashour,  More,  and  Chiu 
[2])  . Stinson,  Davis  and  Khumawala 's  algorithm  performed  well 
for  moderately  sized  job  shop  problems. 
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Talbot  and  Patterson  [77]  described  another  integer 
programming  algorithm  (similar  to  that  of  Balas  [4])  with 
"network  cuts"  to  solve  the  general  resource-constrained 
scheduling  problem.  Instead  of  using  binary  variables,  Talbot 
and  Patterson  represented  the  problem  by  a set  of  structured 
and  compact  arrays.  Their  solution  procedure,  beginning  with 
the  first  of  N jobs,  assigns  to  each  job  its  earliest  feasible 
completion  time  until  either  job  N is  assigned  a completion 
time  or  some  job  j*  (where  j*  < N)  cannot  receive  an 
assignment  due  to  resource  infeasibility.  In  the  latter  case 
backtracking  occurs.  In  particular,  an  effort  is  made  to 
assign  job  j*-l  a completion  time  greater  than  that  previously 
attached  to  it.  In  order  to  fathom  the  partial  solutions 
which  cannot  lead  to  improved  solutions,  an  artifice  called  a 
network  cut  is  introduced.  The  network  cut  is  actually  an 
integer  time  period  which  is  used  to  determine  whether  or  not 
the  associated  partial  schedule  should  be  fathomed.  Talbot 
and  Patterson  showed  that  the  average  CPU  time  reguired  to 
solve  each  of  50  test  problems  from  the  literature  was  reduced 
by  the  employment  of  network  cuts. 

By  the  late  1970s,  there  were  already  more  than  one 
hundred  heuristic  procedures  and  many  exact  solution 
technigues  available  for  RCPSP.  This  motivated  many 
researchers  to  compare  the  solution  technigues  in  an  effort  to 
identify  an  effective  set. 
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In  1975,  Davis  and  Patterson  [13]  compared  the  eight  most 
effective  heuristic  rules  with  an  optimal  solution  procedure, 
developed  by  Davis  and  Heidorn  [12],  called  the  bounded 
enumeration  technique.  The  comparisons  were  performed  on 
eighty-three  test  problems  with  certain  problem 
characteristics.  The  test  results  were  presented  in  the  form 
of  a frequency  distribution  of  the  percent  deviation  from  the 
optimum  for  each  heuristic  rule.  The  minimum  slack  rule 
performed  closest  to  the  optimal  solution  technique.  Finally, 
Davis  and  Patterson  discussed  the  effect  of  problem 
characteristics  upon  the  performance  of  the  heuristic 
procedures.  Patterson  [56]  further  emphasized  the  importance 
of  analyzing  a given  problem's  structure  before  choosing  an 
appropriate  heuristic  to  solve  the  problem. 

Yet  another  experimental  investigation  of  the  performance 
of  heuristics  for  RCPSP  was  conducted  by  Cooper  [9].  He 
examined  the  performance  of  two  heuristic  procedures,  the 
parallel  method  and  the  sampling  method.  The  parallel  method 
generated  a schedule  of  activities  from  a series  of  partial 
schedules  where  a job  was  chosen  to  be  included  according  to 
a particular  priority  rule.  The  sampling  method  selected  the 
best  schedule  among  feasible  alternatives  using  probabilistic 
techniques  that  also  considered  various  job  priority  rules. 
The  projects  that  Cooper  analyzed  in  the  experiments  differed 
in  terms  of  their  size,  structure,  and  resource  properties. 
In  computational  testing  with  the  parallel  method,  he 
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demonstrated  the  superiority  of  the  priority  rule  which  takes 
into  account  all  constraints  when  selecting  the  next  job  to 
commence.  Tests  with  the  sampling  method  failed  to  single  out 
a priority  rule  that  dominates  all  others.  Again,  the  impact 
of  problem  structure  on  the  results  was  noted. 

Patterson  [57]  evaluated  the  three  enumeration-based 
exact  approaches  for  solving  multiple  resource-constrained 
scheduling  problems:  Davis  and  Heidorn's  [12]  bounded 
enumeration  technique,  the  branch  and  bound  technique  of 
Stinson  et  al.  [73]  and  Talbot's  [76]  implicit  enumeration 
technique.  These  techniques  were  compared  using  110  test 
problems  with  a limit  of  five  minutes  of  CPU  time  per  problem 
on  an  AMDAHL  470/V8  computer.  Only  Stinson's  algorithm  solved 
all  110  problems  within  the  specified  CPU  time.  The  effect  of 
problem  characteristics  on  the  solutions  were  examined  by  a 
stepwise  regression  analysis.  Results  demonstrated  the 
superiority  of  Talbot's  algorithm  in  terms  of  computer  storage 
requirements,  Stinson's  algorithm  in  terms  of  computation  time 
needed,  and  Davis  and  Heidorn's  algorithm  in  terms  of 
computation  time  required  when  the  predicted  number  of 
precedence  feasible  subsets  remaining  after  elimination  is 
low. 

A recent  experimental  investigation  was  completed  in  1988 
by  Dumond  and  Mabert  [ 17  ] . They  studied  RCPSP  in  an 
environment  where  new  projects  arrive  continuously  or  randomly 
to  a system  in  which  projects  share  common  resources  and 
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receive  completion  deadlines.  Dumond  and  Mabert  tested  the 
performance  of  four  due  date  procedures  and  five  scheduling 
heuristics  with  full  control  on  the  due  date  assignment.  A 
second  test  was  conducted  to  examine  the  performance  of  the 
due  date  procedures  when  deadlines  were  set  externally.  Their 
experimental  results  failed  to  indicate  a rule  that  uniformly 
outperformed  the  others. 

Levy  and  Tayi  [47]  introduced  the  project  scheduling 
problem  in  a client  contractor  environment.  They  emphasized 
the  cost  overruns  that  occur  as  a result  of  unexpected  project 
redefinitions  during  the  execution  of  the  project.  They 
formulated  a model  that  considered  (a)  the  early, 
intermediate,  and  late  start  options;  (b)  the  project 
crashing  strategies;  and  (c)  the  penalty  costs  associated  with 
project  delays.  They  argued  that  at  certain  stages,  completed 
work  is  evaluated  by  the  client  and  either  accepted  or 
rejected.  It  was  shown  that  the  acceptance  probability  is  a 
nondecreasing  function  of  time  and  has  a large  impact  on  the 
contractor's  optimal  strategy. 

D.  Smith-Daniels  and  V.  Smith-Daniels  incorporated  the 
material  ordering  problem  into  the  project  scheduling  problem 
[72].  Their  formulation  seeks  to  minimize  the  sum  of  the 
costs  due  to  delays  in  completing  the  project,  the  cost  of 
ordering  materials,  the  cost  of  carrying  inventory,  and 
activity  holding  costs.  They  decomposed  the  problem  and 
employed  the  Wagner-Whitin  lot-sizing  algorithm. 
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More  recently,  researchers  have  begun  using  artificial 
intelligence  techniques  in  approaching  project  scheduling 
problems.  Bell  and  Park  [6]  introduced  the  idea  of  using  A* 
search  technique  in  solving  RCPSP.  Partovi  [55]  proposed  a 
set  of  knowledge  based  methods  for  the  selection  and 
application  of  project  management  models. 

Other  recent  research  in  the  area  includes  Weiss  [81], 
Kurtulus  [45]  and  Kurtulus  and  Narula  [46].  Weiss  studied 
RCPSP  in  parallel  networks.  Specifically,  he  formulated  the 
problem  as  an  integer  program  and  solved  the  linear 
programming  (LP)  relaxation  of  his  model  with  the  Dantzig- 
Wolfe  decomposition  technique.  Kurtulus  further  investigated 
the  multi-project  scheduling  problems.  In  particular,  he 
introduced  four  new  scheduling  rules:  (a)  maximum  duration 
and  penalty,  (b)  maximum  penalty,  (c)  maximum  total  duration 
penalty,  and  (d)  maximum  total  work  content.  He  assessed  the 
performance  of  these  as  well  as  six  other  scheduling  rules 
with  respect  to  such  project  summary  measures  as  resource 
constrainedness,  location  of  peak  requirements,  and  problem 
size.  Upon  performing  tests  with  3000  multi-project 
scheduling  problems  with  unequal  and  equal  penalties,  he 
concluded  that  the  maximum  total  work  content  rule  performed 
well  for  small  problems  with  equal  penalties.  Moreover,  he 
found  that  the  maximum  penalty  rule  worked  well  in  solving 
problems  with  unequal  penalties  and  more  constraining  values 
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of  the  average  utilization  factor.  Finally,  it  was  shown  that 
in  all  other  cases  the  minimum  slack  rule  was  most  effective. 

Christofides,  Alvarez-Valdes , and  Tamarit  [8]  suggested 
a branch  and  bound  solution  technique  that  employs  disjunctive 
arcs  to  resolve  resource  conflicts.  Although  they  claimed 
that  this  algorithm  finds  an  optimal  solution  to  a given 
problem,  it  was  later  proven  [15]  with  a counterexample  that 
their  branch  and  bound  procedure  does  not  always  produce  an 
optimal  solution.  However,  was  this  algorithm  modified  to 
obtain  an  optimal  solution  for  the  RCPSP  [14].  Moreover, 
Christophides,  Alvarez-Valdes,  and  Tamarit  added  valid  cuts  to 
the  LP  relaxation  of  the  RCPSP  in  order  to  obtain  strong  lower 
bounds  on  the  project  completion  time  and  compared  these 
bounds  with  optimal  solutions.  They  also  performed 
comparisons  involving  the  Lagrangian  relaxation  of  the  problem 
in  which  the  resource  constraints  were  relaxed.  They  showed 
that  the  Lagrangean  relaxation  yielded  lower  bounds  inferior 
to  those  associated  with  the  LP  relaxation  strengthened  with 
valid  cuts.  Finally,  they  presented  an  alternative  method  for 
obtaining  a lower  bound  for  the  RCPSP.  In  particular,  for  two 
activities  x and  y which  cannot  be  completed  as  a result  of 
resource  conflicts,  they  defined  one  problem  in  which  x 
precedes  y and  another  in  which  y precedes  x.  It  was  shown 
that  one  of  these  problems  provides  a lower  bound  for  the 
original  problem's  optimal  solution. 
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Bell  and  Han  [5],  in  another  recent  study,  described  a 
heuristic  solution  method  for  the  RCPSP.  Their  approach 
consists  of  two  phases;  the  first  step  involves  the  generation 
of  a feasible  schedule,  and  the  second  employs  a "hill- 
climbing" procedure  to  find  local  improvements  in  the  schedule 
found  earlier.  Unlike  other  heuristics,  Bell  and  Han's 
procedure  requires  no  backtracking  in  constructing  a feasible 
solution.  In  particular,  it  proceeds  by  identifying  a 
conflict  or  resource-violating  set  and  then  determining  which 
ordered  pair  of  activities  in  that  set  should  be  linked  by  a 
new  procedure.  Later,  in  the  second  phase  of  the  algorithm, 
the  effect  of  removing  each  of  the  added  arcs  on  the  critical 
path  is  examined  to  improve  the  current  schedule.  In  an 
effort  to  compare  the  performance  of  their  method  versus  that 
of  each  of  several  other  known  heuristics,  Bell  and  Han  coded 
their  procedure  in  Common  Lisp  and  solved  the  110  test 
problems  developed  by  Patterson  [57].  They  discovered  that 
even  though  the  minimum  slack  rule  performed  best  in  terms  of 
the  number  of  problems  solved  optimally  and  the  average 
quality  of  the  near-optimal  solutions,  their  algorithm 
achieved  superior  makespan  results. 

Demeulemeester  and  Herroelen  [14]  have  completed  the 
latest  work  in  the  resource-constrained  project  scheduling 
area.  Specifically,  they  have  suggested  a branch  and  bound 
procedure  which  is  similar  to  a scheme  proposed  by  Johnson, 
[38],  Stinson  et  al.  [73]  and  Christophides  et  al.  [8].  In 
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Demeulemeester  and  Herroelen's  method,  nodes  in  the  depth- 
first  search  tree  correspond  to  feasible  partial  schedules. 
Each  branch  from  the  tree ' s parent  node  corresponds  to  an 
exhaustive  and  minimal  combination  of  activities.  The  delay 
of  these  activities  resolves  a resource  conflict  occurring  at 
the  parent  node.  Two  dominance  rules  are  employed  to  prune 
the  search  tree.  The  first  rule  is  called  the  left-shift 
dominance  rule  and  is  conceptually  identical  to  the  one 
employed  by  Schrage  [66]  and  Stinson  et  al . [73].  The  second 
rule  uses  the  concept  of  a cutset  which  is  defined  as  the  set 
of  all  unscheduled  activities  for  which  predecessor  activities 
are  completed  or  are  still  in  progress. 

Demeulemeester  and  Herroelen  programmed  their  procedure 
in  Turbo  C Version  2.0  on  an  IBM  PS/2  personal  computer. 
Their  method  solved  the  110  Patterson  problems  within  an 
average  of  0.206  seconds.  Furthermore,  their  approach  proved 
to  be  about  12  times  faster  than  the  procedure  of  Stinson  et 
al.  which  was  previously  reported  to  be  the  fastest  for 
solving  the  110  Patterson  test  problems.  In  short,  their 
results  contain  the  best  computational  times  for  minimizing 
the  makespan  of  resource-constrained  project  scheduling 
problems.  For  additional  research  on  this  problem  see  the 
list  of  references  at  the  end  of  this  paper  [66,75]. 

2.2  Time/Cost  Tradeoff  Problem 

As  Phillips  and  Dessouky  [60]  reported,  the  question  of 
primary  concern  in  this  class  of  problems  is  "Which  activities 
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in  a project  should  be  modified  in  duration  in  order  to 
shorten  the  total  project  duration  by  a specified  amount  at  a 
minimum  cost  when  the  duration  of  the  project  exceeds  its 
predetermined  limit?" 

In  1971,  Siemens  [67]  developed  a simple  CPM  time/cost 
tradeoff  algorithm  (SAM) . His  algorithm  starts  with  the 
determination  of  all  paths  from  the  initial  node  to  the  final 
node  and  their  expected  completion  times.  Given  a desired 
project  duration,  his  approach  then  determines  the  amount  by 
which  each  alternative  path  must  be  shortened.  The  effective 
cost  slopes  are  determined  by  modifying  the  actual  cost 
slopes.  The  path  with  the  highest  reduction  and  within  this 
path  the  activity  with  the  lowest  effective  cost  slope  are 
chosen  for  reduction.  Unfortunately,  enumerating  all  paths  in 
the  network  is  computationally  unattractive  even  for 
moderately  sized  problems. 

Phillips  and  Dessouky  [60]  introduced  another  solution 
technique  for  TCTP  by  using  the  minimal-cut  concept.  In  order 
to  achieve  a reduction  in  the  project  duration  at  minimum 
cost,  they  locate  the  minimum  cutset  by  the  out-of-kilter  or 
the  cut  search  algorithm.  They  noted  that  their  algorithm  is 
simple  and  that  problems  can  be  solved  manually.  Furthermore, 
they  remarked  that  the  computational  efficiency  depends  on  the 
cut  algorithm  employed. 

Later,  Goyal  [32]  modified  Siemens'  algorithm  slightly  to 
allow  for  the  deshortening  of  an  activity  duration  that  was 
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previously  shortened.  He  remarked  that  the  problem  Siemens 
faced,  excessively  shortened  paths,  appeared  less  in  the 
revised  algorithm. 

Tufekci  [79]  presented  an  iterative  solution  procedure 
that  employed  a labelling  algorithm  for  locating  a minimum  cut 
in  the  network  flow  representation  of  the  problem  described  by 
Phillips  and  Dessouky  [60].  Through  computational 
experimentation  Tufekci  demonstrated  the  superiority  of  his 
algorithm  over  that  of  Phillips  and  Dessouky. 

Kanda  and  Rao  [41]  incorporated  linear  penalty  costs  for 
delaying  certain  events  into  the  classical  time/cost  tradeoff 
formulation.  Specifically,  they  transformed  the  problem  into 
a sequence  of  algebraic  flow  problems.  The  development  of 
their  proposed  algorithm  is  similar  to  that  of  Fulkerson  [24] 
except  for  modifications  in  the  optimality  conditions  and  the 
introduction  of  the  concepts  of  drained  flow  and  auxiliary 
sources  for  penalty  nodes.  They  tested  their  algorithm  on 
small-scale  problems. 

Hamacher  and  Tufekci' s [34]  work  is  similar  to  that  of 
Kanda  and  Rao  [41]  in  that  the  problem  is  transformed  into  a 
sequence  of  algebraic  flow  problems.  It  is  important  to  note 
that  the  transformation  advanced  by  Hamacher  and  Tufekci  has 
the  ability  to  handle  different  objective  functions.  Further 
research  on  this  problem  is  referenced  at  the  end  of  this 
paper  [ 18 ] . 
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2.3  Payment  Scheduling  Problem 

The  next  type  of  project  scheduling  problem  that  we  will 
survey  in  this  paper  is  the  payment  scheduling  problem.  Given 
a project  described  by  a network,  the  payment  scheduling 
problem  involves  identifying  a schedule  that  maximizes  the 
present  value  of  all  transactions. 

Maximizing  the  net  present  value  of  cash  flows  as  an 
objective  in  a project  with  precedence  relations  was  first 
introduced  by  Russell  [64].  His  formulation  required  the 
maximization  of  a nonlinear  objective  function,  which  is 
neither  concave  nor  convex,  over  a set  of  constraints  that 
represents  the  precedence  relations  between  the  project 
activities.  Russell's  approach  consisted  of  iteratively 
obtaining  a linear  approximation  to  the  objective  function  and 
solving  the  dual  of  the  resulting  linear  program.  The  dual 
linear  program  obtained  in  this  manner  is  a maximum  flow 
problem. 

Subsequently,  Grinold  [33]  demonstrated  that  PSP  can  be 
transformed  into  an  equivalent  linear  programming  problem.  In 
particular,  he  showed  that  the  linear  program  turns  out  to  be 
a weighted  distribution  problem.  He  suggested  two  algorithms: 
the  fixed  deadline  algorithm,  for  solving  the  payment 
scheduling  problem  with  a given  project  deadline;  and  the 
parametric  algorithm,  for  solving  the  problem  for  all  possible 
project  deadlines.  Grinold' s transformation  represents  a 
major  contribution  to  the  project  scheduling  literature 


22 


inasmuch  as  it  continues  to  be  used  in  recent  research  efforts 
in  the  area  (e.g.,  see  Erenguc,  Tufekci  and  Zappe  [19]). 

For  additional  applications  and  extensions  of  the  three 
classes  of  project  scheduling  problems,  see  the  list  of 
references  at  the  end  of  this  paper  [1,10,22,23,31,35,40,  50, 
54,59,71,78]  . 

A considerable  part  of  the  recent  research  in  project 
scheduling  has  been  directed  toward  bridging  the  gap  between 
the  three  classes  of  problems  that  have  been  discussed  above. 

2.4  Payment  Scheduling  with  Cash  Flows  and  Resource 

Restrictions 

As  mentioned  earlier,  most  of  the  work  in  project 
scheduling  has  been  devoted  to  scheduling  subject  to  resource 
constraints  with  the  objective  of  minimizing  project 
completion  time.  In  many  real  life  situations,  however, 
project  managers  are  faced  with  the  problem  of  controlling  the 
cash  flows  and  resource  usage  in  projects  with  hundreds  or 
even  thousands  of  activities.  The  first  efforts  to 
incorporate  the  net  present  value  objective  into  a resource- 
constrained  project  scheduling  problem  were  made  by  Doersch 
and  Patterson  [16]  and  by  Bey,  Doersch,  and  Patterson  [7].  In 
[16]  the  authors  introduced  capital  constraints  into  the 
payment  scheduling  problem.  A set  of  relatively  small  test 
problems  were  solved.  The  computation  time  required  for 
moderately  sized  problems  proved  to  be  excessive.  Later,  in 
1986,  Russell  [65]  studied  the  problem  of  maximizing  the 
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project  net  present  value  subject  to  precedence  relations  and 
resource  restrictions.  He  tested  the  performance  of  six 
heuristic  scheduling  rules  on  this  problem.  Three  of  the 
heuristics  are  often  used  in  unconstrained  cash  flow  analysis, 
one  heuristic  is  a random  rule  which  was  used  as  a benchmark, 
and  the  two  remaining  heuristics  are  commonly  used  in 
minimizing  the  project  duration  subject  to  resource 
constraints.  He  employed  real-life  projects  with  more  than 
one  thousand  activities  in  his  tests.  The  results  showed  that 
for  small  problems  nearly  all  of  the  six  heuristics  performed 
identically.  For  large  scale  problems  with  relatively  loose 
resource  constraints,  however,  the  "minimum  slack  rule" 
performed  the  best.  Finally,  in  solving  problems  with  tight 
resource  constraints  the  superior  performance  was  achieved  by 
the  "target  rule,"  which  assigns  scheduling  priority  on  the 
basis  of  the  difference  between  the  maximal  current  earliest 
finish  time  and  the  unconstrained  optimal  finish  time  for  a 
given  activity  i. 

In  1987,  Smith-Daniels  and  Aquilano  [70]  proposed  the 
idea  of  including  the  project  net  present  value  criterion 
within  "late-start  resource  constrained  project  scheduling." 
They  suggested  an  algorithm  to  determine  whether  or  not  a 
late-start  resource-constrained  schedule,  which  was  shown  to 
be  the  best  schedule  for  minimizing  the  project  duration, 
produces  an  increase  in  project  net  present  value.  The  tests 
performed  on  550  problems  showed  that  there  was  a significant 
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difference  between  the  project  net  present  values  resulting 
from  the  late-start  and  the  early-start  schedules  in  favor  of 
late-start  schedules. 

Recently,  Icmeli  and  Erenguc  [37]  presented  a Tabu  Search 
Algorithm  for  RCPSPDC.  They  tested  their  algorithm  on  50  test 
problems  derived  from  the  Patterson  problems.  The  Tabu  Search 
procedure,  in  general,  produced  near-optimal  solutions  with 
reasonable  computational  effort. 

2.5  Resource  Constrained  Time/Cost  Tradeoff  Project 
Scheduling  Problem 

In  1982,  Talbot  [76]  introduced  the  notion  of  combining 
the  resource-constrained  project  scheduling  and  time/cost 
tradeoff  problems  into  one  formulation  where  the  duration  of 
each  activity  might  vary  according  to  the  resource  allocations 
and  where  the  objective  was  to  minimize  the  project  duration. 
He  presented  an  integer  programming  approach  which  was 
sufficiently  general  to  solve  the  problem  with  various 
objectives,  including  minimal  project  duration,  minimal 
project  cost  given  performance  payments  and  penalties,  and 
minimal  consumption  of  a critical  resource.  An  attractive 
property  of  Talbot's  approach  is  that  his  procedure  can  be 
terminated  with  a good  solution  before  actually  reaching 
optimality.  In  short,  his  approach  produces  optimal  solutions 
for  small  problems  and  good  heuristic  solutions  for  large 
problems . 
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2.6  Time/Cost  Tradeoff  Project  Scheduling  Problem 
with  Discounted  Cash  Flows 

Recently,  Erenguc,  Tufekci  and  Zappe  [19]  integrated  PSP 
and  TCTP.  Given  the  cash  flows  throughout  the  life  of  the 
project,  the  objective  of  the  problem  is  to  determine  the 
activity  durations  and  schedule  of  activity  start  times  so 
that  the  net  present  value  of  cash  flows  is  maximized.  The 
formulation  of  the  problem  enabled  them  to  use  the  Generalized 
Benders  Decomposition  technigue  of  Geoffrion  for  solution. 
Their  algorithm  was  tested  on  56  test  problems.  All  of  the 
test  problems  were  solved  to  optimality  with  modest 
computational  effort. 

2.7  Conclusions 

In  this  chapter,  a unifying  survey  of  the  literature  on 
project  scheduling  is  presented.  The  survey  is  limited  to 
include  the  work  done  specifically  in  the  project  scheduling 
area,  although,  there  are  several  techniques  developed  for 
assembly  line  balancing  and  job  shop  scheduling  that  can  be 
applicable  to  project  scheduling  [53].  It  was  noted  that  much 
of  the  work  in  project  scheduling  was  done  in  the  area  of 
resource  constrained  scheduling  problems.  There  is  some 
recent  research  which  integrates  RCPSP  with  either  TCTP  or 
PSP.  In  addition,  there  is  also  some  work  that  combines  PSP 
and  TCTP.  The  time/cost  tradeoff  problem  and  the  payment 
scheduling  problem  are  linear  programs  and  can  therefore  be 
solved  with  relative  ease.  For  the  third  type  of  fundamental 
project  scheduling  problem,  namely  the  resource  constrained 
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project  scheduling  problem,  there  exist  efficient  exact 
techniques.  Most  notable  among  these  is  the  very  efficient 
branch  and  bound  algorithm  of  Demeulmeester  and  Herroelen 
[14].  However  when  two  or  more  of  these  three  fundamental 
problems  are  integrated,  the  resulting  problems  do  not,  in 
general,  preserve  the  structures  present  in  the  original 
problems.  Consequently,  the  exact  algorithms  available  for 
the  three  fundamental  problems  can  not  be  directly  extended  to 
their  generalizations.  Unfortunately,  in  spite  of  their 
practical  relevance,  for  some  of  the  integrated  problems  no 
research  has  been  done  to  date  and  for  the  others  the 
literature  contains  very  little  work. 

Given  the  financial  relevance  of  the  net  present  value 
objective,  we  expect  to  see  a good  portion  of  the  project 
scheduling  literature  developing  around  the  payment  scheduling 
problem.  Integration  of  PSP,  TCTP  and  RCPSP  and  other 
potential  extensions  of  such  an  integration  appear  to  be  a 
direction  that  the  project  scheduling  literature  is  likely  to 
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CHAPTER  3 

THE  RESOURCE  CONSTRAINED  PROJECT  SCHEDULING  PROBLEM  WITH 

DISCOUNTED  CASH  FLOWS 


3 . 1 The  Problem  Formulation 

The  resource  constrained  project  scheduling  problem  with 
discounted  cash  flows  involves  scheduling  the  project 
activities  with  cash  inflows  and  outflows,  in  such  a way  that 
the  resource  and  precedence  constraints  are  met  and  the  net 
present  values  of  the  cash  flows  that  occur  during  the  life  of 
the  project  is  maximized.  The  following  0-1  ILP  formulation  is 
an  adaptation  of  the  one  proposed  by  Doersh  and  Patterson 
[16].  (P) 

Max  Y Y Qi  e"“fc  xit  ~ Y (t-due)  XNtP  e~at 

i = 1 t =due 


subject  to 


Y xit=1 


t=ef 


Vi 


N t+d1~ 1 

Y Y rixxi<!  * R* 


1=1  q=t 
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Y tx*t  + dn  s £ txnt 


t-e„ 


t=e„ 


( m , n)  eH 
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for  all  (i,j)  e H,  activity  i immediately  precedes 
activity  j 

Xit  e {O,  l)  Vi,t 


where 

T : An  upper  bound  on  the  completion  time  of  the  project. 
Due:  Due  date  on  the  completion  time  of  the  project. 

P : Penalty  per  period  of  delaying  the  project  beyond  its  due 
date. 

N : Last  activity  in  the  project. 

q1  : The  compounded  value  of  cash  flows  that  occur  during  the 
implementation  of  activity  i at  its  completion  time. 

«i  - E cije  ■ 

>i 


dj  : Duration  of  activity  i. 

fjj  : cash  flow  for  activity  i in  period  j,  j=l,...,dj  , f.j 
may  be  negative,  zero  or  positive. 
ei  : Earliest  completion  time  for  activity  i. 
lj  : Latest  completion  time  for  activity  i. 

1 if  activity  i is  completed  in  time  period  t; 

Xit  : ' 

0 otherwise, 
a : Discount  rate. 
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rjk  : Number  of  units  of  resource  k required  by  activity  i. 

Rk  : Availability  level  of  resource  k,  k=l,...,K. 

In  order  to  reduce  the  number  of  integer  variables,  we 
can  determine  the  earliest,  e. , and  latest  completion  times, 
1.,  of  the  activities.  The  method  CPM  is  used  to  determine 
the  earliest  completion  times  of  the  activities.  Thus,  we  do 
not  need  to  consider  the  completion  times  of  activities 
earlier  then  their  ef ' s which  are  determined  by  CPM,  since 
those  will  be  infeasible  completion  times  due  to  the 
precedence  relations.  Similarly,  the  latest  project 
completion  time  can  be  used  to  determine  the  latest  completion 
times  of  the  activities.  The  completion  times  of  the 
activities  later  than  their  lj's  are  infeasible  due  to  the 
precedence  relations. 

The  first  set  of  constraints  requires  that  each  activity 
is  completed  exactly  once.  The  second  set  of  constraints 
ensures  that  resources  allocated  to  activities  at  any  time 
during  the  project  duration  do  not  exceed  the  resource 
availability.  The  third  set  of  constraints  represents  the 
precedence  relations  between  the  activities.  Specifically,  an 
activity  can  only  start  after  all  of  the  activities  which  must 
technologically  precede  that  activity  are  completed. 

The  first  term  in  the  objective  function  is  the  sum  of 
the  values  of  cash  flows  discounted  to  the  beginning  of  the 
project  (time  0)  . The  second  term  is  the  discounted  value  of 
the  delay  penalties  imposed  as  a result  of  delaying  the 
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project  completion  beyond  its  due  date.  The  objective  is  to 
maximize  the  discounted  value  of  all  cash  flows  that  occur 
during  the  life  of  the  project.  Due  date  can  be  determined 
exogenously  or  can  be  set  equal  to  sD,  where  D is  a known 
feasible  project  completion  time  and  s is  a constant  greater 
than  1.  Here,  D is  usually  obtained  from  a heuristic 
procedure.  A one-pass  procedure  for  obtaining  a feasible 
solution  to  RCPSPDC  will  be  presented  in  the  following 
section.  The  latest  completion  time  of  a project,  T,  is 
defined  as  the  due  date  plus  the  maximum  allowable  delay  in 
the  project. 

In  the  model,  the  following  assumptions  are  made: 

1. Each  activity  is  assumed  to 

- have  a known,  deterministic,  integer  duration; 

- consume  known,  deterministic  and  constant  level  of 
resources  of  each  kind  during  its  processing  time; 

- have  no  preemption;  that  is,  once  an  activity  starts, 
it  can  not  be  interrupted. 

2.  The  level  of  availability  of  each  resource  type  is  assumed 
to  be  constant  throughout  the  schedule. 

In  the  following  section,  we  will  discuss  a heuristic 
procedure,  Tabu  Search,  which  recently  gained  practical 
success  in  scheduling  and  other  combinatorial  optimization 
problems.  An  adaptation  of  this  procedure  to  RCPSPDC  will 


follow. 
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3. 2. The  Tabu  Search  Procedure  for  RCPSPDC 
3.2.1.  An  Overview  of  Tabu  Search  Technique 

Tabu  search  was  first  introduced  by  Glover  [27,28]  as  an 
intelligent  search  technique  to  overcome  local  optimality  for 
hard  combinatorial  optimization  problems.  The  simple  Tabu 
search  starts  with  a feasible  solution  and  chooses  the  "best 
move"  according  to  an  evaluation  function  while  taking  steps 
to  assure  that  the  method  does  not  re-visit  a solution 
previously  generated.  This  is  accomplished  by  introducing 
the  Tabu  restrictions  on  the  possible  moves  which  discourage 
the  reversal  and  in  some  cases  repetition  of  selected  moves. 
The  Tabu  list  which  contains  these  forbidden  move  attributes 
is  known  as  the  short  term  memory  function.  This  list  is 
typically  a FIFO  list.  The  most  recent  K moves  or  move 
attributes,  where  K is  the  Tabu  list  size,  are  stored  in  this 
list.  When  a new  member  is  added  to  this  list,  the  member 
recorded  K iterations  ago  is  removed.  The  Aspiration  Criteria 
are  used  to  give  flexibility  to  short  term  memory  process. 
The  Tabu  status  of  a move  can  be  overridden  if  it  meets 
aspiration  criteria. 

Intermediate  and  long  term  functions  may  also  be 
incorporated  to  intensify  and  diversify  the  search. 
Intermediate  term  functions  are  incorporated  by  recording  and 
comparing  the  characteristics  of  best  trial  solutions 
generated  so  far.  The  method  then  seeks  to  generate  solutions 
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which  have  the  same  features  by  restricting  or  penalizing  the 
moves  during  a regional  search. 

The  long  term  memory  functions  guide  the  search  to 
regions  which  are  not  explored  yet.  The  purpose  of  invoking 
a long  term  memory  function  to  search  procedure  is  to  move  to 
a new  region  where  better  solutions  may  be  found  and/or 
entrapment  of  a local  optimal  solution  may  be  avoided. 

Let  T be  the  Tabu  list  and  S(x)  be  the  set  which  consists 
of  moves  s that  can  be  applied  to  x.  Then  the  Tabu  search 
procedure  in  its  simplified  form  is  described  by  Glover  [27] 
as  follows: 

Simple  Tabu  Search 

1.  Select  an  initial  x EX  and  let  x*  = x.  Set  the  iteration 
counter  k=0  and  begin  with  T empty. 

2.  If  S(x)-T  is  empty,  go  to  step  4.  Otherwise,  set  k=k+l  and 
select  sk  ES(x)-T  such  that  sk(x)  = OPTIMUM  (s  (x)  : s ES(x)- 
T)  . 

3.  Let  x=  sk(x)  . If  c(x)  < c(x*)  , where  x*  denotes  the  best 
solution  currently  found,  let  x*  = x. 

4 . If  a chosen  number  of  iterations  has  elapsed  or  since  x* 
was  last  improved  or  if  S(x)-T  = 0 upon  reaching  this  step 
directly  from  step  2,  stop. 

Otherwise,  update  T,  and  return  to  Step  2. 

Figure  3 . Tabu  Search  procedure  in  its  elementary  form 

The  application  of  Tabu  Search  technique  to  optimization 
problems  are  discussed  in  Hertz  et  al.  [36].  After  the 
successful  applications  of  this  approach  to  the  Traveling 
Salesman  Problem  [44,49],  Skorin-Kapov  demonstrated  the 
adaptation  of  the  approach  to  the  Quadratic  Assignment  Problem 
[68],  [69].  Kapov  proposed  a Tabu  Search  procedure,  Tabu- 
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Navigation  Algorithm,  as  a solution  technique  for  QAP.  The 
implementation  of  this  algorithm  allows  user  to  interact 
during  the  procedure.  The  Tabu-Navigation  Algorithm  appeared 
to  be  an  efficient  heuristic  algorithm  for  the  Quadratic 
Assignment  Problem.  The  best  solutions  were  obtained  with 
less  CPU  times.  A comparison  of  the  solutions  obtained  by 
Tabu-Navigation  Algorithm  was  done  with  the  solutions  obtained 
by  Simulated  Annealing  Procedure.  It  was  found  that  Tabu 
search  obtained  better  quality  solutions.  In  her  follow-up 
work  [69],  the  Tabu  search  for  QAP  is  refined  by  including 
knowledge  about  permutations  which  yielded  good  objective 
function  values.  This  method  obtains  the  best  known  results 
for  QAP  problems  with  dimensions  ranging  from  n=42  to  n=90. 

Recently,  Woodruff  and  Spearman  [82]  analyzed  the 
performance  of  Tabu  search  and  simulated  annealing  techniques 
when  these  techniques  are  applied  to  a line  search.  He 
concluded  that  the  time  required  for  move  mechanisms  in  both 
techniques  are  similar  and  the  time  to  actually  make  a move 
does  not  have  an  important  impact  on  the  computation  time. 

Other  Tabu  search  applications  include  the  Tabu  search 
method  designed  by  Gendreau,  Salvail,  and  Soiano  [25]  for  the 
Maximum  Clique  Problem,  by  Fiechter  [21]  for  the  Travelling 
Salesman  Problem,  by  Glover  and  Laguna  [30]  for  the  Bandwidth 
Packing  Problem,  by  Taillard  [74]  for  the  Quadratic  Assignment 
Problem,  and  by  Pirkul  and  Roland  [63]  for  the  Graph 
Partitioning  Problem.  In  the  following  section,  a heuristic 
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procedure  for  RCPSPDC  is  presented.  This  procedure  finds  a 
feasible  solution  for  RCPSPDC  which  will  help  determine  the 
due  date  and  the  upper  bound,  T,  on  the  completion  time  of  the 
project  and  will  constitute  the  initial  feasible  solution  for 
the  Tabu  algorithm  [29,42], 

3.2.2  A Heuristic  Procedure  for  RCPSPDC 

The  heuristic  procedure  starts  with  scheduling  the 
initial  project  activity,  usually  a dummy  activity,  to  start 
at  time  zero  and  proceeds  as  follows.  At  the  completion  time 
of  an  activity,  by  taking  into  consideration  only  the 
precedence  relationships,  all  schedulable  activities  are 
determined.  Let  L be  the  set  of  activities  that  contains  the 
schedulable  activities  at  the  completion  time  of  an  activity. 
If  set  L is  empty,  then  the  procedure  moves  to  the  next  period 
where  an  activity  completion  for  an  already  scheduled  activity 
takes  place.  If  L is  not  empty,  then  the  sum  of  the  resource 
requirements  of  all  the  activities  in  L and  of  those 
activities  that  are  already  in  progress  is  computed  for  each 
type  of  resource.  If  this  sum,  for  each  type  of  resource,  is 
less  than  or  equal  to  the  available  amount  of  resource,  then 
all  of  the  activities  in  L are  scheduled  to  start  at  that 
period.  However,  if  the  total  resource  requirements  exceeds 
the  resource  availability  for  at  least  one  type  of  resource, 
then  a subset  of  activities  in  L is  delayed  until  the  resource 
conflict  is  resolved.  Activities  are  delayed  according  to  the 
following  priority  rules: 
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a.  Delay  the  activity  with  the  largest  (smallest)  cash 
outflow  (inflow) . 

b.  If  rule  (a)  results  in  a tie,  then  delay  the  activity 
with  the  shortest  path  from  itself  to  the  end  of  the 
project. 

The  path  lengths  are  determined  by  taking  into 
consideration  only  the  precedence  relationships. 

Priority  rule  (a)  and  if  necessary  (b)  are  applied  until 
the  resource  conflict  is  resolved. 

The  above  procedure  is  repeated  until  the  last  project 
activity,  usually  a dummy  activity,  is  scheduled.  It  must  be 
clear  that  the  heuristic  procedure  will  always  yield  a 
feasible  solution  to  the  RCPSPDC. 

3.2.3.  A Tabu  Search  Procedure  for  RCPSPDC  (Tabu-S) 

In  order  to  describe  the  Tabu  Search  Procedure  the  terms 
"admissible  move",  "Tabu  List"  and  "evaluation  function"  need 
to  be  defined.  A discussion  of  these  terms  follows. 

In  its  basic  form  a move  is  defined  as  a transition  from 
one  trial  solution  x to  another  solution  x' . The  set  of 
admissible  moves  consist  of  those  moves  that  can  be  applied  to 
the  current  trial  solution,  x.  The  Tabu  List  keeps  record  of 
the  moves  or  move  attributes  which  may  lead  to  re-visiting 
previously  generated  solutions.  Those  Tabu  moves  reside  in 
the  list  for  a specified  number  of  iterations  and  then  become 
free  again.  The  best  move  is  chosen  from  the  admissible  move 
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set  according  to  an  evaluation  function.  The  best  move  is  the 
one  that  takes  the  search  to  a solution  with  the  best 
evaluation  function  value. 

Given  a solution,  a move  is  defined  as  completing  an 
activity  one  time  unit  late  or  early  with  respect  to  its 
current  completion  time.  The  completion  time  resulting  from 
a move  must  be  between  the  earliest  and  latest  completion 
times  of  the  activity  under  consideration.  A move  is  called 
admissible  if  it  is  a non-Tabu  move.  Corresponding  to  each 
activity  there  are  at  most  two  possible  moves  and  consequently 
for  the  entire  schedule  there  are  at  most  2N  admissible  moves 
where  N is  the  total  number  of  activities.  Thus,  the  maximum 
size  of  the  admissible  move  set  is  2N.  A move  may  take  the 
search  to  an  infeasible  solution.  If  a schedule  is  infeasible 
a penalty  is  imposed.  Since  the  problem  is  one  of 
maximization,  this  penalty  has  a negative  sign.  The  sum  of 
the  objective  function  value  of  problem  (P)  and  the  penalty 
,p,  incurred  due  to  infeasibility  is  used  as  the  evaluation 
function,  EF,  for  the  algorithm. 

EF=  £ Yj  Vi  e~at  Xit  ~ £ ( t-due)  XNtP  e"“  t + p 

i=  1 t=e1  t=due 


The  penalty,  p,  associated  with  a solution  is  a 
function  of  the  total  units  of  violation  and  it  is  computed  in 
the  following  manner.  Let  B be  the  total  number  of  precedence 
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and  resource  constraints,  let  s;  be  the  value  of  the  slack 
variable  in  the  ith-  constraint  for  a given  solution.  Then  the 
penalty  associated  with  a particular  solution  is  given  as 

B 

p='mE  vi 

i=l 


where  vt-  = max{  0,-si  } and  M is  a sufficiently  large  positive 
constant. 

As  an  example,  consider  a solution,  x,  which  is  not 
feasible  to  a resource  constraint,  i,  and  a precedence 
constraint,  j,  of  problem  (P) . Then, 

Si  £ rlkxiq)  forsome  k,  t 

1*1  q=t 
In  1* 

Sj  = £ txnt  - { £ txmt  + dn } forsome  activities  m,  n 

t=eD 

Thus, 

s.  <=  0 =»  Vj  = -s.  >=0 

Sj  <=0  =*  Vj  = -Sj  >=0  and 

p = “M  (Vj  + vp. 

The  purpose  of  using  this  penalty  in  the  evaluation 
function  is  first,  to  force  the  feasible  solutions  to  be 
selected  as  the  best  solutions  and  second,  if  there  exists  no 
feasible  solution  resulting  from  the  moves  in  the  set  of 
admissible  moves,  then  the  infeasible  solution  with  the  least 
constraint  violation  is  selected  as  the  best  solution.  This 
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may  potentially  redirect  the  search  towards  the  feasible 
region.  The  best  move  is  a move  in  the  admissible  move  set 
which  takes  the  search  to  a solution  with  the  largest 
evaluation  function  value.  The  corresponding  schedule  may  be 
feasible  or  not.  If  the  current  best  move  takes  the  search 
to  a solution  that  has  a better  evaluation  function  value  than 
all  the  solutions  generated  in  the  previous  iterations,  then 
this  solution  is  kept  as  the  BEST  solution  found  so  far. 

It  is  worth  noting  that  for  a feasible  solution  the 
objective  function  and  the  evaluation  function  values  are  the 
same.  They  differ  only  for  infeasible  solutions. 

Following  is  a formal  statement  of  the  Tabu  Search 
Procedure  for  RCPSPDC: 

1.  Using  the  heuristic  procedure  obtain  a feasible  starting 
solution.  From  this  solution  determine  the  earliest  and 
latest  activity  completion  times. 

2.  Determine  the  Tabu  List  size,  TL  and  the  number  of 
iterations,  MAX. 

3.  At  each  iteration  K (K=l, . . . ,MAX) , 

- Determine  all  the  admissible  moves  from  the  current 
solution. 

- Determine  the  evaluation  function  values  for  all  the 
solutions  resulting  from  admissible  moves. 

- Perform  the  best  admissible  move. 

- Add  the  solution  resulting  from  the  best  move  to  the  Tabu 


List. 
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- If  the  number  of  entries  in  the  Tabu  List  is  greater  than 
TL,  delete  the  oldest  entry  from  the  list. 

- Update  the  best  solution  found  so  far,  BEST,  if 
necessary. 

4.  If  MAX  iterations  have  been  performed,  stop.  Report  BEST 
as  the  best  solution  found. 

It  is  noted  that  if  the  project  due  date  and  the  latest 
completion  time  are  not  exogenously  available,  these  are 
obtained  in  step  1 above  as  discussed  earlier. 

As  it  must  be  clear  from  the  formal  statement  of  the 
Tabu-S  procedure,  rather  than  recording  a solution  itself  in 
the  Tabu  list,  we  chose  to  record  a particular  attribute  of 
it,  its  evaluation  function  value.  Consequently,  in  order  to 
determine  whether  or  not  a particular  solution  is  Tabu,  we 
compare  the  evaluation  function  value  of  the  solution  to  the 
elements  of  the  Tabu  list,  resulting  in  at  most  TL  comparisons 
for  each  potential  move.  However,  had  we  taken  the  approach 
of  recording  the  values  of  the  xjt  variables  in  the  Tabu  list 
and  comparing  the  values  of  the  xit  variables  of  a particular 
solution  to  those  of  the  elements  of  the  Tabu  list,  this  would 
have  required  up  to  (J) (TL)  comparisons  for  each  potential 
move.  Here  J is  the  total  number  of  xjt  variables  in  Problem 
(P)  . Obviously,  our  approach  results  in  a good  deal  of  saving 
in  computational  effort.  This  reduction  in  computational 
effort  comes  at  little  or  no  cost  since  due  to  the  nature  of 
the  evaluation  function  (objective  function) , it  is  unlikely 
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to  obtain  two  or  more  solutions  that  have  exactly  the  same 
function  value. 

The  Tabu  Search  Procedure  for  RCPSPDC  described 
above  makes  use  of  only  the  short  term  memory  function  and  is 
therefore  referred  to  as  Tabu-S.  In  the  following  section  we 
describe  a modification  of  Tabu-S  which  makes  use  of  a long 
term  memory  function  (Glover  [27])  and  is  therefore  referred 
to  as  Tabu-L. 

3.2.4  Introduction  of  Long  Term  Memory  Function  to  Tabu 

Algorithm  (Tabu-L) : 

Preliminary  computational  testing  with  the  Tabu-S 
procedure  indicated  that  the  best  solutions,  BEST,  were 
usually  found  at  the  early  iterations  of  the  algorithm.  This 
motivated  us  to  modify  the  Tabu-S  procedure  to  include  a long 
term  memory  function. 

Tabu-L  is  conducted  in  L stages.  Each  stage  is  started 
with  a different  starting  solution  and  consists  of  m 
iterations.  To  obtain  the  Tabu-L  procedure,  the  following 
modifications  were  made  in  the  Tabu-S  procedure. 

A second  Tabu  list  L-list,  which  records  the  starting 
solution  values  of  each  stage,  is  introduced.  The  purpose  of 
this  list  is  to  prevent  the  algorithm  from  revisiting  any  of 
the  L starting  solutions.  At  each  stage,  the  evaluation 
function  values  of  the  starting  solutions  of  the  current  stage 
and  of  the  previous  stages  are  elements  of  the  L-list.  For 
the  initial  stage,  the  evaluation  function  value  of  the 
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solution  generated  by  the  heuristic  procedure  is  placed  in  the 
L-list.  For  each  stage  1,  1=2,3,...,L,  the  feasible  solution 
generated  in  stage  1-1  that  is  "most  distant"  from  the  best 
solution  obtained  in  stage  1-1  is  selected  as  the  starting 
solution  and  its  evaluation  function  value  is  placed  in  the  L- 
list.  For  any  two  solutions  x and  y the  Euclidian  distance 
between  them  is  used  as  the  distance  measure. 

Let  txi  be  the  completion  time  of  an  activity  i in 
solution  x1  and  ty).  be  the  completion  time  of  the  same  activity 
in  solution  y1 . Then  the  euclidian  distance  between  these  two 
solutions  is  calculated  as 


In  addition,  the  evaluation  function  values  of  the  best 
solutions  obtained  in  each  of  the  L stages  are  made  permanent 
elements  of  the  Tabu  list.  This  prevents  the  procedure  from 
revisiting  any  of  the  best  solutions  obtained  in  the  previous 
stages.  Therefore,  at  any  stage  h,  to  check  whether  or  not  a 
move  is  admissible,  the  evaluation  function  value  resulting 
from  the  move  is  compared  to  the  elements  of  the  Tabu  list  (of 
which  h-1  are  permanent  elements) . By  introducing  the  L-list, 
choosing  the  most  distant  solution  to  the  best  solution  at  a 
stage  as  the  starting  solution  of  the  next  stage,  and  making 
the  best  solution  of  each  stage  a permanent  member  of  the  Tabu 
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list  we  aim  to  diversify  the  search  procedure  to  regions  that 
markedly  contrast  with  those  examined  so  far  [27]. 

Following  is  a formal  statement  of  the  Tabu-L  procedure: 

1.  Using  the  heuristic  procedure  obtain  a feasible  starting 
solution  for  the  first  stage  (1=1) . From  this  solution 
determine  the  earliest  and  the  latest  activity  completion 
times.  Add  the  objective  value  of  the  solution  to  the  L-list. 

2.  Determine  the  Tabu  list  size,  TL,  the  number  of  iterations 
at  each  stage,  m and  the  number  of  stages,  L. 

3.  At  each  stage  1 (1=1,...,L): 

3.1.  At  each  iteration  K,  (K=l,...,m): 

3.1.1.  Determine  all  the  admissible  moves  from  the 
current  solution. 

3.1.2.  Determine  the  evaluation  function  values  for  all 
the  solutions  resulting  from  admissible  moves. 

3.1.3.  Perform  the  best  admissible  move. 

3.1.4.  Add  the  solution  value  resulting  from  the  best 
move  to  the  Tabu  list. 

3.1.5.  If  the  number  of  entries  in  the  Tabu  list  is 
greater  than  TL,  delete  the  oldest  entry  from 
the  list. 

3.1.6.  Update  the  best  solution  found  so  far  if 
necessary . 

3.2.  At  the  completion  of  m iterations,  define  as  the 
"most  distant"  solution,  the  most  recently  generated 
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feasible  solution  whose  objective  value  is  not  a 
member  of  the  L-list. 

3.3.  Add  the  best  solution  value  of  stage  1 to  Tabu  list 
as  a permanent  member. 

3.4.  Let  the  "most  distant"  solution  be  the  starting 
solution  of  stage  1+1.  Add  this  solution  value  to 
L-list. 

4.  If  L stages  have  been  performed  report  the  BEST  solution 
found  so  far. 

Note  that  in  step  3.2  of  each  stage  1,  an  approximation 
to  the  Euclidian  distance  is  used  for  determining  the  "most 
distant"  solution  to  the  best  solution.  If  the  best  solution 
for  a given  stage  is  obtained  in  the  early  iterations,  which 
is  usually  the  case,  then  this  approach  will  yield  a "good" 
approximation . 

3.2.5  Computational  Results 

The  two  Tabu  Search  Procedures  were  coded  in  Fortran 
and  tested  on  50  project  networks  selected  from  the  Patterson 
[17]  data  set.  The  computational  experiments  were  done  on  an 
IBM  3090-600J  Computer  with  vector  processing.  In  order  to 
determine  the  values  of  the  parameters  TL,  MAX  and  m,  several 
preliminary  runs  were  performed.  Increasing  the  Tabu  list 
size  beyond  10  did  not  affect  the  solution  quality.  Therefore 
the  Tabu  list  size  was  set  at  10.  As  indicated  earlier,  in 
the  implementation  of  Tabu-S,  the  best  solutions  were  obtained 
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rather  early  in  the  procedure.  Our  preliminary  runs  indicated 
that,  for  the  problem  dimensions  we  are  dealing  with,  in  a 
majority  of  the  problems  the  best  solution  was  obtained  within 
the  first  hundred  iterations  and  furthermore,  in  no  case  was 
the  best  solution  obtained  later  than  the  370th  iteration. 
These  observations  and  the  fact  that  we  wanted  to  make  the 
experiments  with  Tabu-S  and  Tabu-L  comparable  provided  the 
justification  for  letting  L=4  and  m=100  and  MAX=400.  Since  in 
the  Tabu-L  procedure  the  best  solutions  obtained  at  each  stage 
are  permanent  members  of  the  Tabu  list,  this  list  may  contain 
up  to  14  moves.  (Note  that  at  the  first  stage  Tabu  list  has 
at  most  11  members,  at  the  second  stage  it  has  at  most  12 
members , etc . ) 

Project  due  date  and  the  upper  bound,  T,  on  the 
completion  time  of  the  project  are  computed  as  follows: 

Due  date=  project  completion  time  found  by  the  heuristic 
procedure* 1 . 25 
T = Due  date  + s 
where 


2,  if  due  date  <=10 


s= 


3,  otherwise 

Since  the  due  dates  on  the  completion  times  of  the 
projects  are  not  readily  available,  we  approximated  the  due 
dates  by  using  the  project  completion  times  found  by  the 
heuristic  procedure  plus  a tolerance  to  make  the  problems  as 
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realistic  as  possible.  To  determine  a feasible  upper  bound  on 
the  project  completion  time,  a tolerance,  s,  was  added  to  the 
due  date.  To  make  the  upper  bound  more  meaningful,  tolerance 
for  longer  projects  was  made  slightly  larger. 

The  activity  cash  flows  were  generated  randomly  from  a 
uniform  distribution  on  [-500,1000].  The  penalties  were 
assigned  as  follows: 

r500  per  time  period,  if  due  date  < 10 

P = , 

1000,  otherwise 

The  purpose  of  making  P larger  for  longer  projects  is  to  take 
into  account  that  discounting  will  make  the  cash  outflows, 
that  are  further  in  the  future,  less  penalizing. 

There  are  no  special  purpose  exact  solution  techniques 
for  the  RCPSPDC  to  the  authors'  knowledge.  Attempts  were  made 
to  solve  several  of  the  test  problems  using  a general  purpose 
ILP  software  (LINDO) . Problems  with  7 and  9 activities  were 
solved  for  an  optimal  solution  on  the  IBM  3090-600J  mainframe 
within  approximately  60  CPU  seconds.  However,  for  problems 
with  22  activities,  a feasible  solution  could  not  be  obtained 
within  200  CPU  seconds.  Consequently,  we  did  not  experiment 
with  LINDO  for  the  larger  problems.  Since  the  optimal 
solutions  for  these  problems  are  not  readily  available,  the 
solutions  obtained  by  Tabu-S  and  Tabu-L  were  compared  with 
upper  bounds  obtained  from  the  linear  programming  relaxations 
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of  RCPSPDC  that  are  strengthened  with  valid  inequalities. 

3. 2. 5.1.  An  Upper  Bounding  Scheme  for  RCPSPDC 

An  upper  bound  on  the  value  of  any  integer  (or  mixed 
integer)  Linear  Program  can  be  determined  by  solving  the 
Linear  Programming  relaxation  of  that  problem.  Thus,  for 
maximization  problems,  the  optimal  solution  to  any  integer  (or 
mixed  integer)  linear  programming  yields  a value  less  than  or 
equal  to  the  value  of  the  optimal  solution  obtained  by 
relaxing  its  integrality  requirements. 

If  the  integrality  requirement  of  problem  (P)  was 
relaxed,  then  the  following  formulation  is  obtained: 

(P* ) 

Max  E E e~Bt  xit  - E It-due)  XNtP  e~at 

i=  1 t=ei  t=due 

subject  to 

E ^t=1 

t=ei 
jy  t+dj-1 

E E ri*xi«  * R* 

i=l  g=C 

i.  In 

E tx*t  + dn  i E ^ 

t=eD 

for  all  (i,j)  e H,  activity  i immediately  precedes 


Vi 


Vic,  fc 


activity  j 
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Xit  ;>  0 Vi,  t 

(P1)  is  a linear  problem  which  can  be  solved  with  any 
commercial  software  package.  The  value  of  the  optimum  solution 
obtained  by  solving  problem  ( P * ) gives  an  upper  bound  on  the 
optimal  value  of  the  problem  (P) . However,  the  gap  between 
these  two  solutions  may  be  large.  In  order  to  strengthen  the 
linear  programming  relaxation,  ( P • ) , secondary  constraints  or 
valid  inequalities  can  be  added.  These  inequalities 

effectively  cut  (or  eliminate)  certain  parts  of  the  solution 
space  that  do  not  contain  feasible  integer  solutions. 

Consequently,  the  optimality  gap  can  be  reduced. 

In  general,  a valid  inequality  can  be  defined  as 
follows. 

Given  a polyhedron  S={x  e Rn  : Ax  <=  b },  The  inequality 
nx  <=  nQ  is  called  a valid  inequality  for  S if  it  is  satisfied 
by  all  points  in  S.  Thus,  irx  <=  nQ  is  a valid  inequality  if 

and  only  if  S lies  in  the  half  space  {x  e Rn  : nx  <=  nQ  } or 

equivalently  if  and  only  if  max  { nx:  x e S } <=  nQ. 

Proposition  3.1  : (I)  is  a sets  of  valid  inequalities  for 

problem  (P) , where 

Sy  ={ 1 , . . . , i, . . . , j , . . . ,N)  and  S ={Sy  : v=l,...,V),  i.e.  S is 
the  set  of  all  paths  in  the  project. 
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t+c^-l  t+dj-1  t+d„-l 

E ^iU+---+  E ^+---+  E xJq+---+  E 

u=t  m=t  g=t  l=ew 

Vt,  VSv  e5  ( eg.  I) 

Proposition  3.2  : (II)  is  a set  of  valid  inequalities  for 

problem  (P)  where 

t+dj-l 

E E xiv  * min*{L*}  Vt  {eg.  II) 

ielt  q=t 


and  where  It  contains  the  activities  that  can  be  active  at 
time  period  t. 

The  first  set  of  inequalities  (I)  is  based  on  the 
precedence  relations  between  the  activities  on  a path  of  a 
project.  A path  is  defined  as  a sequence  of  arcs  (or 
activities)  in  which  the  initial  node  of  each  arc  is  the  same 
as  the  terminal  node  of  the  preceding  arc  in  the  sequence 
where  nodes  correspond  to  the  beginning  and  ending  of  the 
activities.  Thus,  by  definition,  an  activity  on  a path  can 
not  be  done  simultaneously  with  any  other  activity  on  the  same 
path.  This  conclusion  led  us  develop  the  first  set  of  valid 
inequalities  (I) . 

Let  Sy  denote  a path  in  a network  where  l,i,j  and  N are 
the  activities  on  this  path.  Let  S denote  the  set  of  all 
paths  in  the  project.  If  activity  1 completes  at  time  t then 
this  indicates  that  it  is  active  at  time  t.  But  other 
activities , i , j,  N,  on  the  path  can  not  be  completed  at  time 


49 


t through  time  t plus  their  durations,  since  that  indicates 
that  they  are  active  at  time  t which  is  not  possible  by  the 
above  path  definition.  If  activity  1 is  the  initial  activity 
of  the  path,  then  all  the  activities  on  the  path  must  succeed 
this  activity.  This  further  implies  that  if  an  activity  is 
active  at  time  t,  then  activity  1 must  be  completed  before 
time  t and  it  can  not  be  active  at  any  time  grater  than  or 
equal  to  time  t.  Thus,  the  first  summation  in  the  inequality 
may  go  up  to  1,.  Similarly,  if  N is  the  last  activity  on  the 
path,  then  all  the  activities  on  the  path  must  precede  this 
activity.  Thus,  activity  N can  not  be  scheduled  to  complete 
at  its  earliest  completion  time  unless  all  other  preceding 
activities  on  the  path  are  scheduled  to  their  earliest 
completion  time.  Therefore,  the  last  summation  may  start  from 

eN- 

The  second  set  of  constraints  (II)  is  based  on  resource 
utilization.  Let  It  be  the  set  of  schedulable  activities  at 
time  t which  can  not  be  done  simultaneously,  since  scheduling 
all  the  activities  in  It  causes  resource  conflict  for  some 
resource  type  k.  Lk  denotes  the  maximum  number  of  activities 
that  can  be  scheduled  without  the  resource  violation  at  that 
time  period.  Then  (II)  implies  that  Lk  or  less  activities 
which  are  in  It  can  only  be  scheduled  to  complete  at  any  time 
between  t and  t plus  their  durations,  since  completing  an 
activity  at  that  time  period  denotes  that  activity  is  active 
at  time  t.  Since  for  each  resource  type  k the  set  It  will  be 
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the  same  for  some  time  t,  the  inequality  with  mink  (Lk  } is 
sufficient  for  consideration. 

A pseudo-algorithm  for  determining  the  set  It  and  Lu: 

Do  1 for  each  time  t=l,...,T 

It  = 0 

Do  2 for  each  activity  i=l,...,N 
If  ej-dj+1  <=  t <=  lj,  then  It  <-  It  u {i} 

2 Continue 

1 Continue  { set  It  for  all  t has  been  found  } 

Do  3 for  resource  type  k=l,...,K 

Sort  the  activities  in  It  in  an  increasing  order  of  their 
resource  requirements,  r(X).  It  ={ X(1) , X(2) , . . . , X(M) } where 

r ( x<  i > ) <=  r(x(2>)  <=.••<=  r(X(M)). 

Sum=0 

Count=l 

Do  4 for  i=l , . . . , N 
Sum  = Sum  + r(X(i)) 

If  Sum  > Rk 

Then  Lk=  Count- 1 go  to  3 
Else  count  = count+1 
4 Continue 

3 Continue  { Lk  is  found  } 

Given  the  above  explanations,  the  validity  of  these  two 
inequalities  can  be  proved. 
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Proof  of  Proposition  3.1  : (I) 

Let  X be  the  feasible  region  defined  by  constraint  sets 
(1) , (2) , (3)  and  (4)  of  problem  (P) . Furthermore,  let  x be 

a solution  of  (P)  such  that  xu1  = x2t2  = ...=  xNtM  = 1.  Let 
1,  i,  j , k,  1 , . . . ,p,N  be  activities  on  a path  Sy  of  (P)  . Then  by 
definition  of  a path,  1 « i « j « k « 1 «...«  p <<  N, 
where  k <<  1 means  activity  k immediately  precedes  activity  1. 
Since  x e X,  by  (3)  , the  following  relations  hold  for  the 
activities  on  a path  Sy:  (Rl) 

t.  - tl  >=  d.  =>  t,  >=  d.+t, 


*"N  ^ > dN  > ^"N  > dN+tp 

Then,  by  above  relations,  if  activity  i completes  at  time 
t,- , then  t1  <=  t.-d.  and  tN  >=tp>=.  . . >=tj>=d.+tj . This  implies 
that  at  time  tj-d^l  through  t.+dj-l,  none  of  the  activities  on 
the  path  can  be  completed.  Thus,  at  time  t1 , 

T xiu=  E • • • = E E *m  = o, 

u=CJ  g=  Cj  v=tt 

ci+di_1 

E ^=1 


Since  activity  i and  path  Sy  are  chosen  arbitrarily  xeX 
must  satisfy  the  following  inequality. 
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1'en 

(eg.  I) 


Vfc,  V5v  eS 


Q . E . D. 


Proof  of  Proposition  3.2:  (II) 

Let  X be  the  feasible  region  defined  by  the  constraint 
sets  (1)  , (2)  , (3)  and  (4)  of  problem  (P)  . Let 


where  It  contains  the  activities  that  can  be  active  at  time 
period  t. 

Show  that  x e X =>  x e X2. 

Let  x be  a solution  of  X where  x1tl  = x2t2  =...=  xNtN  = 1. 
Thus,  x satisfies 


That  is,  solution  x is  a resource  feasible  solution.  Let  Zt, 
be  the  set  of  activities  which  are  active  at  time  t'  in  this 
solution  for  some  resource  type  k.  Thus, 


t+^-i 


*2  = { X X xiq  * min^L*}  Vt,lc  } 


ielt  q=t 


If  t+c^-1 


X X rikXiq  * Rk 


i=  1 q=  C 
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t'+dj-l 

E E rikXiq*  Rk  <*> 

ieZt/  q=t> 

Let  the  number  of  activities  in  Z{1,  be  Lk'.  By  the 
construction  of  Lk,  given  in  the  above  algorithm,  we  know  that 
at  most  Lk  activities  can  be  active  without  resource  violation 
in  any  feasible  solution  at  time  t'.  Since  Zt,  is  resource 
feasible  by  (k)  , we  must  have  Lk  >=  Lk'  or 

t'+dj-l 

E E xi<,  * v s Lk 

ieZc/  gr=  t1 


As  explained  before,  for  each  resource  type  k,  the  set  Zt,  will 
be  the  same  at  time  t',  thus,  it  will  be  sufficient  to 
consider  the  inequality  with  mink  { Lk  } , since  the  other  valid 
inequalities  will  be  redundant  at  time  t'.  But  t'  is  chosen 
arbitrarily.  Therefore,  x e X2.  Q.E.D. 

The  total  number  of  valid  inequalities  derived  above  is 
T * (total  number  of  paths  in  a network)  + T.  Thus  for  large 
problems,  the  number  of  constraints  may  be  substantial.  In 
order  to  reduce  the  dimension  of  the  constraint  set,  a subset 
of  the  set  of  paths  in  the  network  can  be  considered. 

3. 2. 5. 2 Conclusions 

The  above  upper  bounding  scheme  is  used  to  obtain  upper 
bounds  for  the  test  problems.  In  driving  valid  inequalities 
we  chose  the  paths  according  to  the  following  rule: 
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- Starting  with  the  path  containing  the  maximum  number  of 
activities,  we  selected  paths  until  each  activity  of  the 
network  was  included  in  at  least  one  path. 

The  added  inequalities,  in  general,  reduced  the  upper 
bounds  obtained  from  the  LP  relaxations  of  the  problems  by  1 
to  5%. 

Computational  results  are  summarized  in  table  1.  In 
columns  3 through  7,  the  first  and  second  entries  in  each  row 
correspond  to  results  with  0.01  and  0.02  discount  rates, 
respectively.  (It  is  noted  that  in  his  computational 
experiments  Russell  [64]  used  a discount  rate  of  0.01.) 

In  terms  of  the  gap  between  the  solution  obtained  and  the 
upper  bound,  Tabu-L  performed  at  least  as  well  as  Tabu-S  in  41 
of  the  50  problems.  In  the  remaining  9,  Tabu-S  did  slightly 
better.  In  terms  of  computational  efforts  Tabu-S  required 
slightly  less  CPU  time.  Accordingly,  the  following 
discussions  will  be  limited  to  Tabu-L.  Following  are  the 
average  CPU  times  consumed  by  Tabu-L  for  the  problems 

solved. 

Table  1.  The  average  CPU  times  consumed  by  the  problems  with 
the  indicated  number  of  activities. 

number  of  activities  CPU  times  (seconds) 

10  to  26  6.0 

27  8.8 

35  17.4 
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46.0 
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We  recognize  the  fact  that  these  CPU  times  may  be 
improved  with  more  effective  coding. 

For  the  50  problems  solved,  Tabu-L  produced  solutions 
that  were,  on  the  average,  within  0.03  and  0.05  of  optimality 
for  0.01  and  0.02  discount  rates  respectively.  The  ranges  of 
optimality  gaps  were  [0.001,0.200]  and  [0.002,0.270]  for 
a=0.01  and  a=0.02,  respectively.  Ninety  eight  percent  of  the 
problem  with  a=0.01  and  eighty  percent  of  the  problems  with 
a=0.02  yielded  solutions  that  were  within  ten  percent  of 
optimality.  It  must  be  noted  that  the  optimality  gap  analysis 
discussed  above  is  based  on  an  upper  bound  not  on  the  optimal 
value  of  the  problem.  Accordingly,  one  would  expect  the 
optimality  gap  to  be  smaller.  It  is  also  noted  that  the 
optimality  gap  is  larger  for  problems  with  a=0.02.  This  is 
not  necessarily  an  unfavorable  indication  on  the  part  of  the 
Tabu  search  procedures,  it  may  be  at  least  partly  due  to  the 
weakening  of  the  linear  programming  relaxations. 

In  the  literature,  several  studies  demonstrated  that  the 
Minimum  Slack  (MINSLK)  heuristic  performed  well  for  RCPSP 
[13,57].  In  resolving  resource  conflicts,  the  Minimum  Slack 
Heuristic  gives  priority  to  the  activity  with  minimum  slack 
which  is  determined  by  the  critical  path  analysis  [65]. 
Russell  [65]  computationally  tested  the  performance  of  several 
heuristics  for  RCPSPDC.  He  recommended  the  MINSLK  heuristic 
for  small  to  medium  size,  moderately  resource  constrained 
problems.  Accordingly,  MINSLK  solutions  for  the  50  test 
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problems  were  obtained  and  included  in  the  computational 
table.  Tabu-L  obtained  better  solutions  in  39  of  the  50  test 
problems  with  a=0.01  and  in  42  of  the  50  test  problems  with 
a=0.02  (see  Table  2).  As  one  would  expect  Minimum  Slack 
heuristic  performed  slightly  better  for  small  problems  with 
activities  most  of  which  have  cash  inflows. 


TABLE  2.  Comparison  of  Tabu-L,  Tabu-S  and  Minimum  Slack 

Heuristic  solutions  with  the  solutions  obtained  from 
the  LP  relaxation  of  RCPSPDC  strengthened  with  the 
strong  valid  inequalities . 


Problem 

Number  of 
Activities 

Minimum 

Slack 

Heuristic 

Tabu-S 

solution 

Tabu-L 

solution 

LP  with 
valid 
ineq. 

% Gap 

1 

14 

30171.77 

30150.17 

30217.53 

30314.89 

0.003 

27030.33 

27108.76 

27215.37 

27269.80 

0.002 

7 

17150.91 

17144.34 

17144.34 

17185.62 

0.002 

2 

16644.03 

16543.28 

16543.28 

16676.36 

0.008 

13 

20284.80 

21154.65 

21154.65 

21732.52 

0.03 

3 

17135.17 

18083.75 

18083.75 

21162.12 

0.14 

22 

23179.70 

23187.70 

23195.62 

23284.71 

0.004 

4 

22310.20 

22325.04 

22353.85 

22498.06 

0.006 

22 

1068.86 

1093.37 

1114.66 

1157.45 

0.002 

5 

1038.05 

1084.80 

1100.11 

1210.02 

0.09 

22 

4094.32 

4114.77 

4114.77 

4126.64 

0.003 

6 

3901.68 

3935.16 

3940.85 

3962.91 

0.005 

9 

22225.91 

22204.00 

22212.07 

22285.55 

0.003 

7 

21015.65 

20972.83 

20972.83 

21158.26 

0.008 

9 

7290.46 

7463.72 

7463.72 

7746.32 

0.03 

8 

6410.39 

6725.71 

6727.34 

7036.42 

0.04 

18 

6790.21 

6552.12 

6556.93 

7287.55 

0.10 

9 

5873.76 

5470.14 

5525.46 

6685.25 

0.17 

10 

10071.00 

10093.12 

10096.80 

10114.69 

0.001 

10 

9063.85 

9102.14 

9108.37 

9137.77 

0.003 

11 

8 

356.84 

356.84 

356.84 

374.33 

0.04 

306.56 

306.56 

306.56 

335.38 

0.08 

Table  2 — continued 
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Problem 

Number  of 
Activities 

Minimum 

Slack 

Heuristic 

Tabu-S 

solution 

Tabu-L 

solution 

LP  with 
valid 
ineq. 

% Gap 

12 

23 

4874.68 

5005.66 

5005.66 

5037.48 

0.006 

4366.06 

4610.89 

4610.89 

4666.66 

0.01 

13 

22 

2228.67 

2160.73 

2127.58 

2568.17 

0.17 

1911.10 

1780.18 

1770.51 

2440.12 

0.27 

14 

35 

6454.20 

6578.40 

6595.10 

7320.56 

0.11 

5206.65 

5419.10 

5435.60 

6142.22 

0.13 

15 

35 

1641.85 

1708.07 

1708.07 

1969.74 

0.13 

847.70 

1027.00 

1026.13 

1334.82 

0.23 

16 

22 

2314.69 

2472.06 

2472.09 

2669.74 

0.05 

1958.77 

2112.26 

2209.00 

2515.96 

0.12 

17 

22 

1877.57 

2116.05 

2116.05 

2186.88 

0.03 

1325.81 

1417.72 

1417.72 

1753.50 

0.19 

18 

22 

2663.06 

2816.05 

2802.82 

2975.33 

0.05 

2220.92 

2417.72 

2422.72 

2714.89 

0.10 

19 

22 

3474.73 

3613.96 

3590.10 

3689.18 

0.02 

2745.82 

2978.07 

2978.07 

3084.46 

0.03 

20 

22 

1596.66 

1645.25 

1644.36 

1696.59 

0.03 

1245.32 

1328.26 

1325.19 

1404.28 

0.05 

21 

22 

3511.58 

3580.44 

3569.12 

3911.76 

0.08 

2590.01 

2699.43 

2672.03 

3072.44 

0.13 

22 

22 

3495.86 

3677.30 

3677.30 

3892.04 

0.05 

2709.36 

2977.84 

3012.08 

3344.78 

0.09 

23 

22 

5169.86 

5366.80 

5366.80 

5732.58 

0.06 

3845.55 

4140.16 

4140.16 

4716.67 

0.12 

24 

22 

4379.68 

4502.91 

4519.21 

4654.78 

0.03 

3252.82 

3464.94 

3481.57 

3620.83 

0.04 

25 

22 

5625.10 

5577.41 

5580.66 

5986.13 

0.06 

4590.48 

4526.64 

4531.15 

5170.18 

0.12 

26 

22 

4802.91 

4842.09 

4844.28 

4935.37 

0.01 

4196.00 

4189.13 

4258.79 

4420.08 

0.03 

27 

22 

3893.38 

4047.09 

4048.83 

4290.88 

0.06 

3159.60 

3388.51 

3388.51 

3659.59 

0.08 

28 

22 

5113.66 

5381.23 

5381.23 

5473.85 

0.01 

4317.54 

4763.06 

4763.06 

4918.99 

0.03 

Table  2 — continued 
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Problem 

Number  of 
Activities 

Minimum 

Slack 

Heuristic 

Tabu-S 

solution 

Tabu-L 

solution 

LP  with 
valid 
ineq. 

% Gap 

29 

22 

3864.08 

3832.64 

3838.11 

3959.58 

0.03 

3406.57 

3353.70 

3361.78 

3583.57 

0.06 

30 

22 

3165.09 

3186.14 

3190.00 

3427.13 

0.06 

2563.46 

2691.34 

2690.59 

3013.71 

0.10 

31 

27 

5972.35 

6472.87 

6472.87 

6769.30 

0.04 

4459.96 

5130.75 

5148.50 

5653.85 

0.08 

32 

27 

6387.04 

6597.35 

6614.07 

6780.14 

0.02 

5114.03 

5434.16 

5434.16 

5718.10 

0.05 

33 

27 

8819.46 

8791.09 

8790.04 

9096.19 

0.03 

7589.49 

7646.44 

7646.44 

8037.46 

0.04 

34 

27 

8572.48 

8744.26 

8744.26 

9040.44 

0.003 

7166.79 

7447.28 

7447.28 

7668.30 

0.03 

35 

27 

4693.47 

4677.25 

4615.20 

5116.10 

0.10 

3735.99 

3614.29 

3612.06 

4045.44 

0.12 

36 

27 

5717.42 

5902.98 

5934.06 

6171.42 

0.03 

4634.92 

4633.35 

4930.48 

5177.00 

0.05 

37 

27 

5429.10 

5483.84 

5483.84 

5537.83 

0.01 

4392.33 

4428.12 

4439.60 

4567.75 

0.03 

38 

27 

5925.05 

6159.51 

6170.30 

6355.10 

0.03 

4646.48 

5349.10 

5359.10 

5734.20 

0.07 

39 

27 

4839.23 

5010.16 

5041.32 

5391.61 

0.06 

3847.65 

4150.35 

4153.06 

4435.20 

0.06 

40 

27 

5655.19 

5894.48 

5888.46 

6477.30 

0.10 

4775.98 

5033.45 

5130.55 

5694.30 

0.11 

41 

51 

9097.48 

9022.23 

9063.29 

9878.67 

0.09 

6692.38 

8675.12 

8251.99 

9571.16 

0.16 

42 

22 

4558.77 

4722.55 

4722.55 

4958.65 

0.03 

3889.79 

4168.03 

4168.03 

4501.44 

0.08 

43 

22 

6096.40 

6109.97 

6109.97 

6218.56 

0.02 

5270.92 

5010.24 

5010.24 

5490.80 

0.08 

44 

22 

4489.16 

4514.00 

4539.85 

5446.80 

0.20 

3632.17 

3912.98 

3712.98 

4565.76 

0.23 

45 

22 

4606.80 

4647.30 

4622.46 

4899.23 

0.06 

4046.16 

4107.66 

4106.07 

4598.72 

0.13 

Table  2 — continued 
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Problem 

Number  of 
Activities 

Minimum 

Slack 

Heuristic 

Tabu-S 

solution 

Tabu-L 

solution 

LP  with 
valid 
ineq. 

% Gap 

46 

51 

9566.84 

10274.65 

10277.17 

10379.77 

0.01 

7050.33 

9034.23 

9466.17 

%55.49 

0.02 

47 

51 

11036.71 

6932.12 

9738.14 

10614.42 

0.09 

8213.21 

8831.72 

8889.60 

9778.56 

0.10 

48 

51 

12070.11 

9265.40 

9265.40 

12324.52 

0.05 

9464.15 

8441.23 

8441.23 

9837.% 

0.09 

49 

51 

8684.86 

10542.62 

10542.62 

11174.52 

0.06 

6102.30 

9718.88 

9719.06 

104%.58 

0.08 

50 

27 

5795.29 

5894.68 

5902.64 

6085.44 

0.001 

4852.21 

5077.% 

5076.59 

5228.80 

0.03 

% Gap-(LP  with  valid  ineq.-  Tabu_L  sol.)/LP  with  valid  ineq. 


3.2.6  Concluding  Remarks 

In  this  chapter,  two  Tabu  search  procedures  for  the 
resource  constrained  project  scheduling  problem  with 
discounted  cash  flows  are  proposed.  Tabu-L  uses  a long  term 
memory  function  and  Tabu-S  relies  on  a short  term  memory 
function.  Computational  testing  indicated  that  of  the  two 
approaches,  Tabu-L  would  be  the  preferred  one  for  the  problems 
solved.  Both  procedures  were  tested  on  50  test  problems. 
Computational  results  demonstrated  Tabu  search  as  a viable 
approach  for  the  RCPSPDC  in  that  for  all  the  problems  solved, 
"good  quality"  solutions  were  obtained  with  reasonable 
computational  effort. 

In  the  following  section,  an  exact  solution  procedure  is 


suggested  for  RCPSPDC. 
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3.3  A Branch  and  Bound  Procedure  for  RCPSPDC 
3.3.1  Introduction 

In  this  section,  a branch  and  bound  procedure  is 
developed  for  RCPSPDC.  Although  there  are  many  exact  solution 
techniques  for  RCPSP,  to  our  knowledge  there  does  not  exist 
any  exact  solution  techniques  for  RCPSPDC.  We  are  motivated 
by  the  opportunity  for  providing  a set  of  benchmark  problems 
with  their  optimal  solutions  which  can  be  used  for  testing  the 
quality  of  exact  and  heuristic  solution  procedures  to  be 
developed  in  the  future. 

The  branch  and  bound  technique  is  the  most  prominent 
search  method  for  integer  programming  problems.  It  originates 
from  the  straightforward  idea  of  enumerating  all  feasible 
integer  points.  The  basic  idea  is  to  develop  clever  tests 
that  consider  only  a (small)  portion  of  the  feasible  region 
explicitly  but  automatically  account  for  the  remaining  points 
implicitly. 

The  following  results  will  be  needed  for  the  development 
of  the  branch  and  bound  procedure  for  RCPSPDC. 

Let  (P1')  be  the  problem  obtained  by  removing  the 
resource  constraints  from  the  constraint  set  of  problem  (P) 
and  by  including  the  penalty  value  of  objective  function  as  a 
cash  outflow,  of  the  last  activity,  of  a network. 

We  will  obtain  an  optimal  solution  to  problem  ( P ' ' ) by 
solving  a payment  scheduling  problem.  Therefore,  to  be 
consistent  with  the  PSP  literature  we  will  change  our  notation 
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used  in  problem  (P)  and  problem  (P1')  slightly.  Let  each 
activity  be  denoted  by  (i,j)eA  where  there  are  N activities  in 
set  A.  Let  the  last  activity  of  the  project  be  (s,L). 
Furthermore,  let  V={1,...,L}  correspond  to  events,  Tf  be  the 
realization  time  of  event  i and  e;j.  and  1^.  be  the  earliest  and 
latest  completion  times  of  the  activity  (i,j). 

Let  qn-  j ' denote  the  modified  cash  flows  on  the 
activities.  That  is, 

f qf j,  if  (i/j)  * (s,L), 


l>  j 


.<3jj  + P(t-due)  , if  (i,j)  = (s,L)  and  t > due 


Then,  (P ' ' ) 


1u 

Max  E E <*i}‘  e'at  xijt 

(i.j)eA  t*eL] 


subject  to 


■Lij 

E xut=1  V (i,  j)  eA 


(eq.l) 


tl 


— YU 

E tXmvt  + dvu  * E tXvut  ( (m,v)  , (V,U))  EH 


t=e. 


t=e„ 


(eq. 2 ) 


ijt 


e (o,  l)  V (i,  j)  eA,  V t 


(eq.3) 
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The  payment  scheduling  problem  can  be  formulated  as  follows: 
(PSP) 

M 

Max  f(T)  ='jT/  cie~aTl 

i= 1 

subject  to 

Tj  -T±  S:  d±j  for  j =2 , 3 , . . ,M  , iePj 
T1Z  0 i =1,  ..  ,M 

where 

Tj  : Realization  time  of  event  i,  i=l,..,M.  T^O. 

c.  : Cash  flow  at  the  realization  of  event  i. 

djj  : Duration  of  activity  (i,j)  e A ( A is  the  set  of 

activities  in  the  project.)-  d,-j's  are  integer  valued. 
a : Discount  rate. 

Pj  : Set  of  events  immediately  preceding  event  i. 

Fi  : Set  of  events  immediately  following  event  i. 

(j  is  in  Pj  if  and  only  if  (j,i)  e A;  j is  in  Fj  if  and 
only  if  (i,j)  e A.) 

The  activity  (M,l)  links  the  completion  event  M and 
initial  event  1.  DH1  is  an  upper  bound  on  the  completion  time 
of  the  project,  i.e.,  DHl  = T in  problem  (P)  . 

From  Grinold  (33],  given  integer  djj's,  event  occurrence 
times,  Tj , i=l, . . . ,M,  are  integer  valued. 

Proposition  3.3.1:  Problem  ( P ' ' ) can  be  written  as  a 

payment  scheduling  problem  (PSP) . 
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Proof ; Given  the  latest  completion  time  of  the  project,  the 
earliest,  e{ , and  the  latest  occurrence  times,  1.,  of  each 
event  can  be  determined.  The  earliest  and  the  latest 
occurrence  times  of  events  can  be  determined  by  the  critical 
path  method. 


Let  X 


it 


'1,  if  event  i is  realized  at  time  t 


.0,  otherwise 


Then  PSP  can  be  written  as:  (PSP') 

Max  E E cis"Cxit 

i- 1 t-ei 

subject  to 

U ly 

Y txuc  + duv^Y  txvt  V ( u,  vj  eA. 

t=ev 

Y xit  = 1 Vi,i  = i,...,M. 

t=ei 

Xlce  { 0 , l } Vi,  t 


(eq. 2a) 


(eq.2b) 
(eq. 2c) 


Problem  (P' ')  involves  determining  the  completion  time  of 
each  activity  such  that  the  net  present  value  of  cash  flows 
which  occur  at  the  completion  times  of  the  activities  is 
maximized  subject  to  precedence  constraints  and  subject  to  the 
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constraints  which  guarantee  that  each  activity  completes 
exactly  once. 

Problem  (PSP1)  involves  scheduling  the  occurrence  time  of 
each  event  such  that  the  net  present  value  of  cash  flows  on 
the  events  is  maximized  subject  to  precedence  constraints  and 
subject  to  the  constraints  which  guarantee  that  each  event 
occurs  exactly  once. 

Claim:  Let  (i,j)  be  an  activity  with  > 0 (q^.  < 0).  Then 
at  an  optimum  solution  of  problem  (P' ') , activity  (i,  j) 
starts  (completes)  at  Tf  (Tj)  . 

Proof  of  claim:  Easily  follows  from  the  fact  that  when  q^. 

> 0,  then  qi- j e'at  > q^.  e'at‘  where  t < t'.  Similarly,  when  q^. 

< 0,  then  q(j.  e'at  < q^.  e'at‘  where  t > t'.  Q.E.D. 

Then,  by  the  above  claim,  we  can  relate  the  objective 
function  of  problem  ( P * ' ) to  that  of  (PSP1): 

ci=  ^2  qije~adii  where  qij*  0 

(i.j)  eFj 

+ ^2  Qij  where  q-j  < 0 . (z) 

(.i.j)  ePt 

The  above  relation  shows  that,  cash  flows  on  the  arcs 
(activities)  can  be  equivalently  represented  by  cash  flows  on 
the  nodes  (events) . 

Note  that  each  event  represents  the  completion  of  all 
activities  leading  into  it.  Therefore,  the  precedence 
relations  between  the  activities  can  be  precisely  expressed  in 
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terms  of  the  precedence  relationships  between  the  events, 
thus,  constraints  (2)  can  be  written  as  constraints  (2a)  . 
Constraints  (1)  along  with  constraints  (3)  assumes  that  each 
activity  completes  exactly  once  and  by  definition  of  an 
"event"  this  can  be  equivalently  stated  as  each  event  must 
occur  exactly  once.  Hence,  constraints  (1)  and  (3)  can  be 
written  as  constraints  (2b)  and  (2c).  Q.E.D. 

Upon  solving  PSP,  the  optimal  activity  start  times  (and 
completion  times)  can  easily  be  obtained  from  the  optimal 
event  occurrence  times,  ' s.  That  is,  by  the  above  claim, 
if  activity  (i,j)  has  non  negative  cash  flow,  then  at  the 
optimal  solution  it  will  be  scheduled  to  start  at  time  T-*. 
Similarly,  if  activity  (i,j)  has  negative  cash  flow,  then  it 
will  complete  at  time  Tj* . 

In  his  paper,  Grinold  [33]  has  shown  that  the  payment 
scheduling  problem  is  equivalent  to  a linear  programming 
problem  and  that  the  problem  has  an  optimal  solution 
corresponding  to  a spanning  tree,  an  extreme  point  in  the  set 
of  feasible  schedules, in  the  project  network.  This  special 
structure  is  used  to  develop  two  algorithms  where  the  first 
one,  the  fixed  deadline  algorithm,  solves  the  problem  for  a 
given  deadline  and  the  second  algorithm  finds  optimal 
solutions  for  all  possible  deadlines  . The  fixed  deadline 
algorithm  is  explained  in  the  appendix  A.  This  efficient 
solution  procedure  is  used  in  our  branch  and  bound  procedure 
to  solve  the  subproblems  at  each  node  optimally. 
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3.3.2  The  Search  Process 

At  the  initial  node  of  the  tree  an  optimal  solution  to 
the  (P'1)  problem  is  obtained  by  solving  the  corresponding 
(PSP)  problem.  The  objective  function  value,  UB°,  is  an  upper 
bound  for  the  (P)  problem.  If  this  solution  is  resource 
feasible,  then  it  is  also  an  optimal  solution  for  the  problem 
(P)  and  the  procedure  terminates  at  this  level.  If,  on  the 
other  hand,  there  is  at  least  one  period  with  resource 
conflict,  then  we  branch  from  this  solution  to  resolve  that 
conflict.  The  resource  conflicts  are  resolved  by  introducing 
additional  precedence  relations  to  problem  (P'').  The 
subproblems  thus  obtained  are  solved  as  payment  scheduling 
problems.  Associated  with  each  node  there  is  a set  of 
additional  precedence  relations  imposed  on  the  problem  (P11). 
Let  Ck(p)  be  the  set  that  contains  all  the  activities  which 
are  active  at  the  time  when  the  conflict  occurs  at  node  k. 
The  delaying  set  Dk(p)  at  level  p consists  of  all  subsets  of 
activities  Dkq,  q=l,...,Q,  the  delay  of  which  would  resolve 
that  particular  conflict  within  the  activities  in  Ck(p) . The 
delaying  alternative  Dkq  is  minimal  if  it  does  not  contain 
other  delaying  alternatives  as  a subset  (similar  definitions 
used  in  Demeulemeester  et  al.  [14]).  Each  delaying 
alternative  Dq  is  scheduled  at  the  completion  time  of  some 
activity  in  Ck(p)  . The  delay  of  Dkq  is  introduced  by  defining 
additional  precedence  relations  between  the  activities  of  Dkq 
and  Ck(p)\Dkq. 
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In  order  to  facilitate  the  presentation  of  the  branch  and 
bound  algorithm  we  will  introduce  the  following  notation. 

t+djj-l 

H:{xijt'  £ £ rijk  Xijq±Rk  Vk.t) 

(i,j)6A  q=t 

1 mv  1 vu 

GiiXij  t:  £ tXmvt  + d^,  s Y,  txvut  Um.v)  ,{v,u))eH  .xljt*  0,V. 

t~emv  t=evu 

1u 

I : bcijt : Xijt  = 1 , e {o , l)  , V ( i , j)  EA,  V t } 

t=eij 

Thus,  we  can  write  problem  (P)  as  max  f(x) 

s.t.  x e GnHnI 

and  problem  ( P 1 • ) as  max  f(x) 

s.t.  x e GnI. 

Bpk  : Total  number  of  branches  generated  at  node  k 
at  level  p of  the  tree. 

Gk  : The  set  of  precedence  relations  that  the  network 
has  at  the  kth-  node  of  the  tree.  These 
include  the  original  precedence  and  those  that  have 
been  added.  (Gk=G°  = G and  the  nodes  of  the  tree 
are  numbered  in  the  order  they  are  generated.) 

GPNk  : The  set  of  precedence  relations  at  the  parent  node 

of  node  k. 

LB  : The  best  solution  found  so  far  for  problem  (P)  in 
the  branch  and  bound  tree,i.e.,  the  lower  bound  for 
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(P) . LB  is  set  equal  to  the  value  of  a feasible 
solution  if  one  is  available. 

NB  : The  set  of  nodes  which  are  not  branched  from  yet. 
P' ' (k)  : The  subproblem  solved  at  each  node  k of  the  tree. 

UBk  : Optimal  objective  function  value  found  at  node 

k. 

Assume  that  we  are  at  node  k of  level  p of  the  branch  and 
bound  tree.  We  obtain  an  optimal  solution  to  problem  P' 1 (k) , 

P' ' (k)  max  f (x) 

s.t.  x e Gknl 

by  solving  the  corresponding  PSP.  The  optimal  solution  to 
P1  ' (k)  , UBk' , gives  an  upper  bound  on  the  optimal  solution  UBk 
of  P(k) , 

P(k)  max  f(x) 

s.t.  x e GknHnI. 

Let  xk  denote  the  optimal  solution  to  P' ' (k) . Clearly  xk  e 
Gknl  if  furthermore  xk  e GknHnI. 

Note  also  that  the  resource  constraints  in  problem  (P) 
can  be  implied  by  additional  set  of  precedence  relations,  GA. 
That  is, 

(P)  max  f(x) 

s.t.  xe  Gm  n I 

where  GM  = G°n  GA. 

Let  UB*  be  the  optimal  solution  to  (P) . Then, 

UB°  > UB*  > LB. 


rel .3.3.2 
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As  will  be  discussed  shortly,  at  any  node  k,  the  number 
Bpk,  of  branches  that  the  algorithm  can  generate  is  at  most 
N(N-l).  Each  branch  corresponds  to  delaying  a set  of 
activities  whose  delay  prevents  the  particular  conflict,  which 
initiated  the  branching  operation,  from  occurring  again.  At 
these  newly  generated  nodes  the  precedence  relationships,  to 
accomplish  the  specific  delay  objectives,  are  imposed  and  the 
resulting  problems  P' ' (i) 's  i=l,2, . . . , Bpk  are  solved.  For  each 
problem  solved  a check  is  made  to  see  if  the  node  can  be 
fathomed  and/or  the  incumbent  can  be  updated.  If  a node  can 
not  be  fathomed  then  it  is  added  to  the  candidate  list,  NB. 
The  algorithm  proceeds  by  choosing  a node  (subproblem)  from 
the  candidate  list  and  by  repeating  the  branching  and 
fathoming  operations  discussed  above.  The  algorithm 
terminates  when  the  candidate  list  becomes  empty. 

Proposition  3.3.2  : In  order  to  resolve  a particular  resource 

conflict,  it  is  sufficient  to  consider 
only  minimal  delaying  alternatives. 

Proof : Assume  that  at  a branch  of  some  level  we  have  a 
delaying  alternative  Dq  which  is  not  minimal.  Assume  with  out 
loss  of  generality  that  there  is  an  extra  activity  in  the  set 
Dq.  Let  z'  be  the  optimal  solution  obtained  when  the 
precedence  relations  are  added  to  the  problem  due  to  the  delay 
of  Dq.  Furthermore,  let  z"  be  the  optimal  solution  obtained 
when  this  extra  activity  is  not  delayed.  Assume  to  the 
contrary  that  z'>z".  But  this  is  not  possible  because  if 
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delaying  that  activity  is  profitable,  then  when  we  solve  the 
problem  without  imposing  a precedence  relation,  the  optimal 
solution  delays  that  extra  activity  indifferently. 
Consequently,  z'  <=  z".  Q.E.D. 

More  Definitions: 

f(x)  : The  objective  function  of  problem  (P) . 

N : Total  number  of  activities  in  a project. 

N (A)  : Total  number  of  elements  in  set  A. 

p : A level  of  branch  and  bound  tree. 

PRpb  : The  set  of  precedence  relations  imposed  to 

problem  (P1 ')  at  branch  b,  to  solve  the  resource 
conflict  at  level  p of  the  tree.  b=l,..,Bpk. 

PN  : Parent  node. 

xINC  : Incumbent  solution. 

3.3.3  The  Branch  and  Bound  Algorithm 

Step  0.  Let  p=0 , k=0,  n=0,  LB=-  ,NB=0.  Solve  problem 
P' 1 (0) : max  f(x)  s.t.  x e G°  n I.  If  the  optimal 
solution,  xQ  e GnHnI  then,  stop  and  return  UB°; 
xQ  is  an  optimal  solution  to  problem  (P) ; 
else  continue;  Add  node  n to  the  candidate  list; 
Let  p <-  (p+1) . 

Step  1.  Determine  the  first  resource  conflict  of  the 
schedule  obtained  at  node  n. 

1.1.  If  no  resource  conflict  occurs,  then  go 
to  step  2.1. 
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Else,  determine  Cn(p)  and  Dn(p)  . Branch 
from  node  n to  resolve  this  conflict. 
That  is,  determine  Bpn  where, 

Bp=  £ N (C  n (p)  \Dn2) 

iT\ 

and  PR  b. 

p 

Step  2.  Corresponding  to  each  branch  at  level  p, 

generate  a node  k, where  the  problem  P''(k) 
is  solved:  Max  f(x)  s.t.  x e Gknl,  where, 

Gk=  GPNknPRpb.  Choose  a node  n such  that, 

UBn  = max  { UBk:  k e DPN(p)  }. 

2.1.  Check  all  nodes  and  fathom  those  with 
UBk  < LB.  If  UBn  < LB,  then  fathom  all 
nodes  at  level  p and  go  to  step  3 ; Else 
if  the  solution,  xn  e GnlnH,  and  UBn  >= 
LB  then,  LB  <-  UBn,  xINC=xn.  Fathom  all 
the  nodes  in  DPN(p)  Go  to  step  3;  Else 
xn£GnInH,  add  nodes  in  DPN(p)\n  to  the 
candidate  list,  NB,  p <-  (p+1) , go  to 
step  1. 

Step  3.  Backtrack.  Set  p <-  ( p— 1 ) . 

3.1.  If  p=l,  then,  go  to  step  4.  Else 
continue. 

3.2.  Choose  an  unfathomed  node  at  this  level 
such  that,  UBn  = max  {UBk  : ke  DPN(p)  and 
keNB  }. 
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3.3.  If  there  exist  no  unfathomed  node  at  this 
level,  set  p <-  (p-1)  and  go  to  step  3.1. 
Else,  go  to  step  1. 

Step  4.  Termination.  Return  LB,  which  is  the  optimal 
solution  value  to  problem  (P) . 

Note  that  the  initial  LB  for  the  problem  can  be  obtained 
by  a heuristic  procedure. 

Claim  1:  The  maximum  number  of  branches  generated  at  each  node 
of  the  branch  and  bound  tree  at  each  level  is  N(N-l)  . 
Proof : Each  branch  solves  a resource  conflict  by  adding 

additional  precedence  relations  between  the  activities  at  some 
level  p.  By  the  definition  of  Bpk  ,by  (1),  the  maximum  number 
of  activities  in  set  Cn(p)  is  N.  Furthermore  the  minimum 
number  of  activities  in  the  delaying  alternative  Dnt  is  1. 
Therefore,  N(Dk(p))=  N and  N (Ck  (p)  \Dkq)  = N-l  for  all  q=l,..,N 
and 

N 

Bp=J^  (N-l)  ^N(N-l) 

l=i 


Claim  2 : The  maximum  number  of  precedence  relations  in  a 
project  network  is  N(N-l)/2. 

Proof : Assume  that  we  have  a1#...,aN  activities  in  a project 
network.  Then  activity  a,  can  have  maximum  (N-l)  precedence 
relations  with  the  activities  a2, . . . ,aN.  Since  project  networks 
are  directed  networks,  activity  a2  can  have  maximum  (N-2) 
precedence  relations  when  activity  a1  has  (N-l)  precedence 
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relations.  Therefore,  for  N activities  in  a project  network, 
there  will  be  a total  of  (N-l)  + (N-2) + . . .+1  = N(N-l)/2 

precedence  relations  that  can  be  defined.  Q.E.D. 

Proposition  3.3.3:  The  above  branch  and  bound  algorithm 

finds  the  optimal  solution  to  problem  (P)  , 
in  a finite  number  of  steps. 

Proof ; It  should  be  clear  by  section  3.3.3  that  the  algorithm 
terminates  when  the  candidate  list,  NB,  becomes  empty.  The 
incumbent  solution,  xINC,  is  an  optimal  solution  to  (P)  . This 
follows  from  the  validity  of  upper  and  lower  bounds  given  in 
(rel.3.3.2)  and  it  further  follows  from  the  fact  that  at  each 
level  of  the  branch  and  bound  tree  we  consider  all  possible 
ways  of  resolving  a particular  resource  conflict  by  adding 
precedence  relations  between  the  activities  causing  this 
conflict.  Then,  one  of  the  branches  at  a given  level  should 
contain  the  precedence  relation (s)  that  is  required  by  the 
optimal  solution  that  the  algorithm  will  yield.  Thus,  no 
feasible  solution  to  (P)  is  excluded  from  consideration. 

In  this  second  part  of  the  proof,  we  now  show  that  the 
algorithm  must  terminate  after  a finite  number  of  stages. 

The  algorithm  terminates  when  NB=o.  If  NB  becomes  empty 
before  all  the  terminal  nodes  are  generated  then,  trivially, 
the  algorithm  terminates  in  a finite  number  of  stages. 
Nontrivial  case. 

By  claim  1,  the  maximum  number  of  branches  generated  at 
each  node  of  the  branch  and  bound  tree  at  each  level  is  N(N- 
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1) . Furthermore  by  claim  2,  the  maximum  number  of  precedence 
relations  a project  network  can  have  is  N(N-l)/2.  Thus,  the 
maximum  number  of  nodes  that  can  be  generated  is  [N(N-l)]2  /2. 
Let  K c[N  (N-l)  ]2/2  be  the  number  of  terminal  nodes.  From  the 
argument  given  above,  it  is  clear  that  K is  finite.  Let  NB 
contain  these  K nodes.  Let  x*  be  an  optimal  solution  to  (P)  . 
Then  clearly  there  exist  a terminal  node  n*  whose  optimal 
value  is  f(x*),  optimum  objective  function  value.  Let  na  be 
any  terminal  node  and  xa  be  the  optimal  solution  obtained  at 
this  node.  Then,  f(xa)  < f(x*)  and  thus  na  is  fathomed  from 
the  candidate  list.  Since  na  was  an  arbitrary  element  of  the 
candidate  list,  the  proof  is  complete.  Q.E.D. 

3.3.4.  An  Example 

This  section  gives  an  incomplete  example  of  the  branch 
and  bound  procedure  for  RCPSPDC.  The  first  two  levels  of  the 
tree  will  be  formed  which  will  be  sufficient  for  reader  to 
understand  the  key  concepts  of  the  procedure.  Initially  there 
are  seven  activities  and  6 events.  The  deadline  of  the 
project  is  set  to  12.  The  discount  factor  is  0.01  monthly. 
The  following  table  summarizes  the  data  for  the  problem. 
Assume  that  we  have  one  resource  type  and  Rgv  =3 . 

Table  3 

Activity  cjj  dj  Immed.  Successor 


1 

2 

3 

4 

5 

6 
7 


100  1 

-200  2 

-700  2 

400  3 

700  2 

-200  1 

500  1 


5 
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Note  that  q.  • s are  the  cash  flows  discounted  to  the  nodes. 


Figure  3. Project  Network  and  the  Branch  and  Bound  Tree 

At  k=0  and  p=0,  the  optimum  solution  to  the  relaxed  problem  is 
UB°=  502.37.  Since  this  solution  is  not  resource  feasible 
we  determine  the  first  resource  conflict:  p=l,  C°(l)={2,3}  and 
D°(l)={  {2}, {3}  }.  Thus,  we  have  two  branches  running  from 

node  0 to  solve  resource  conflict:  delay  3 after  2 and  delay 
2 after  3.  Two  more  nodes  are  generated  where  in  the  first 
one  we  solve  the  relaxed  problem  with  the  additional 
precedence  relation  between  2 and  3 and  in  the  second  node  we 
solve  the  problem  with  the  additional  precedence  relation 
between  3 and  2.  In  the  first  node  we  have  obtained  UB1  = 
475.  This  solution  is  not  resource  feasible.  Then,  at  p=2  , 
C1  ( 2 ) ={  3 , 4 , 7 } and  d’(2)=  { {3}, {4}, {7}  }.  The  total  number 

of  branches  generated  from  this  node  will  be  6.  In  the  second 
node,  we  have  obtained  UB2  = 482  and  this  solution  is  not 
resource  feasible  either.  C2(2)={4,7)  and  D2(2)=  {{4}, {7}}. 

There  are  two  branches  running  from  this  solution.  Since  UB2 
> UB1  , we  continue  from  the  second  node.  p=3  and  at  the 
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third  node  we  solve  the  problem  with  an  additional  precedence 
relation  between  4 and  7.  UB3=  4 34.1  and  the  solution  is 
resource  feasible.  Update  the  lower  bound,  LB=434.1.  And  the 
procedure  continues  by  solving  the  problem  at  node  4. 

3.3.5  Computational  Results 

The  Branch  and  Bound  Algorithm  was  coded  in  FORTRAN  and 
tested  on  the  same  set  of  project  networks  that  was  utilized 
for  testing  the  performance  of  the  Tabu  Search  Procedures  (see 
section  3.2.5) . The  computational  experiments  were  done  on  an 
IBM  3090-600J  computer  with  vector  processing.  For  each 
problem  the  maximum  number  of  nodes  generated  was  set  to  2000 
and  the  maximum  computer  processing  time  was  set  to  30 
minutes.  The  feasible  solutions  generated  by  the  heuristic 
procedure  described  in  section  3.2.2  were  used  as  the  initial 
lower  bounds  for  the  problems. 

For  the  problems  with  less  than  51  activities  we  use  a 
tolerance  of  0.5%.  That  is,  a node  with  an  upper  bound  which 
is  within  0.5%  of  the  current  LB  is  fathomed.  Preliminary 
testing  indicated  that  for  the  problems  with  51  activities, 
the  Branch  and  Bound  algorithm  terminated  without  finding  a 
feasible  solution  within  the  given  CPU  limit.  Consequently, 
for  problems  with  51  activities  the  tolerance  was  increased  to 
5%.  Tables  4 and  Table  5 give  a summery  of  computational 
testing  done  on  the  50  test  problems  with  0.01  and  0.02 
discount  rates  respectively.  In  each  table,  the  second  column 
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demonstrates  the  best  solutions  found  with  0.5%  tolerance  for 
the  problems  with  less  than  51  activities  and  with  5% 
tolerance  for  the  problems  with  51  activities.  The  total 
number  of  subproblems  solved  for  each  problem  is  given  in 
column  3.  Column  4 shows  the  total  number  of  simplex  (Grinold 
[33])  pivots  done  for  each  problem.  The  average  number  of 
pivots  for  each  subproblem  is  given  in  column  5.  Note  that 
rather  than  solving  a new  subproblem  from  scratch,  the  Branch 
and  Bound  procedure  updates  the  solution  obtained  at  the 
parent  node  of  a newly  generated  node  to  reoptimize  at  this 
new  node.  Therefore,  it  can  be  claimed  that  columns  3 and  4 
demonstrate  the  computational  efficiency  of  the  problem. 
Column  6 demonstrates  the  maximum  number  of  children  (maximum 
size  of  a conflict  set)  that  any  node  had,  and  column  7 gives 
the  average  number  of  children  of  a tree.  Column  8 gives  the 
total  number  of  precedence  relations  added  at  the  optimal 
solution.  Note  that  these  additional  precedence  relations  are 
required  to  resolve  resource  conflicts.  Column  9 demonstrates 
the  maximum  list  size  which  indicates  the  maximum  number  of 
unfathomed  nodes  at  any  stage  of  the  procedure.  Note  that  the 
maximum  list  size  is  the  size  of  NB. 

Finally,  the  last  column  gives  the  CPU  times  in  seconds. 
The  following  two  tables,  Table  6 and  Table  7,  summarize  the 
results  in  the  preceding  two  tables.  From  Tables  6 and  7,  the 
following  observations  are  evident: 
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- Columns  3,4,5  and  9 clearly  demonstrate  the  efficiency  of 
the  algorithm:  The  average  number  of  nodes  generated  to  find 
an  optimal  solution  in  each  group  is  significantly  lower  than 
the  maximum  limit  of  2000.  The  average  number  of  pivoting 
done  at  each  of  these  nodes  for  each  group  of  problems  is  at 
most  6.  Lastly,  the  average  number  of  unfathomed  nodes  at  any 
point  during  the  procedure  is  at  most  103 . 

- Overall  performance  of  the  algorithm  did  not  change  as  a 
result  of  a change  in  the  discount  rates. 

- As  the  problem  size  increased  (number  of  activities) , the 
number  of  precedence  relations  added  to  obtain  an  optimal 
solution  and  the  size  of  the  branch  and  bound  tree  (depth  and 
width  of  the  tree)  did  not  expand  significantly. 

- For  problems  with  dense  precedence  relations  and/or  with 
large  number  of  activities  CPU  times  increased  moderately. 

It  was  also  noted  that  determining  the  sets,  U and  V of 
the  spanning  tree  in  the  Grinold  (see  appendix  A)  procedure 
used  for  solving  subproblems  accounted  for  a big  portion  of 
the  total  CPU  time.  Therefore,  with  more  effective  coding  and 
improved  way  of  determining  the  sets  U and  V,  these  CPU  times 
can  be  improved  significantly. 

Table  4.  B&B  results  with  0.5%  tolerance  and  discount 


rate  of  0.01 


Problem 

Optimum 

Solution 

Totl.  no 
subprob . 

Total  no 
piv/prob 

Ave  no 
piv/subpr 

max  confl 
set 

Ave.  no 
children 

No  of 
preced 
added/pr 

Max  list 
size 

CPU 

time 

(sec) 

1 

30276.33 

27 

85 

2 

6 

5 

3 

6 

264 

2 

17156.83 

3 

8 

2 

2 

3 

2 

2 

0.42 
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Table  4 — continued 


Problem 

Optimum 

Solution 

Totl.  no 
subprob . 

Total  no 
piv/prob 

Ave  no 
piv/subpr 

max  confl 
set 

Ave.  no 
children 

No  of 
preced 
added/pr 

Max  list 
size 

CPU 

time 

(sec) 

3 

21618.07 

92 

112 

1 

9 

4 

4 

26 

8.97 

4 

23284.71 

1 

8 

8 

0 

0 

0 

1 

0.65 

5 

1157.46 

1 

11 

11 

0 

0 

0 

1 

0.76 

6 

4123.05 

13 

100 

7 

12 

13 

2 

13 

6.43 

7 

2225.20 

15 

34 

1 

6 

5 

2 

10 

1.60 

8 

7631.23 

12 

48 

2 

3 

2 

3 

7 

1.66 

9 

7251.44 

1055 

5638 

4 

24 

8 

5 

63 

265 

10 

10114.66 

1 

3 

3 

0 

0 

0 

1 

0.22 

11 

374.32 

3 

4 

1 

2 

3 

2 

3 

0.41 

12 

4945.18 

1023 

5432 

4 

15 

3 

17 

83 

873 

13 

2544.10 

1314 

9927 

6 

19 

5 

4 

48 

640 

14 

6643.32 

442 

2700 

5 

20 

9 

3 

35 

675 

15 

1966.47 

181 

1546 

7 

12 

6 

3 

31 

387 

16 

2663.97 

453 

3345 

6 

8 

6 

5 

33 

219 

17 

2144.24 

2052 

7374 

3 

12 

7 

11 

71 

582 

18 

2938.58 

115 

820 

5 

6 

3 

5 

25 

55 

19 

3598.96 

106 

306 

2 

12 

6 

17 

90 

398 

20 

1661.59 

73 

293 

3 

12 

5 

10 

49 

593 

21 

3825.20 

1918 

9662 

3 

9 

2 

16 

26 

788 

22 

3892.04 

67 

337 

4 

6 

6 

2 

27 

24 

23 

5517.34 

39 

45 

1 

7 

3 

13 

29 

341 

24 

4557.08 

1129 

2894 

2 

17 

4 

21 

66 

468 

25 

5778.89 

737 

2327 

2 

15 

9 

14 

94 

540 

26 

4895.45 

732 

1931 

2 

12 

7 

6 

47 

467 

27 

4103.51 

1776 

6124 

2 

20 

4 

18 

99 

590 

28 

5449.66 

1859 

4075 

1 

12 

6 

9 

85 

440 

29 

3924.30 

586 

2653 

3 

8 

3 

5 

34 

210 

30 

3280.06 

112 

401 

3 

7 

4 

5 

22 

356 

31 

6676.44 

128 

576 

4 

16 

8 

17 

109 

1191 

32 

6754.05 

67 

64 

1 

8 

6 

12 

57 

337 

33 

8945.77 

1505 

9610 

5 

12 

6 

13 

89 

1432 

34 

8864.53 

241 

621 

2 

12 

6 

13 

80 

777 

80 


Table  4 — continued 


Problem 

Optimum 

Solution 

Totl.  no 
subprob. 

Total  no 
piv/prob 

Ave  no 
piv/subpr 

max  confl 
set 

Ave.  no 
children 

No  of 
preced 
added/pr 

Max  list 
size 

CPU 

time 

(sec) 

35 

4825.92 

125 

512 

3 

12 

7 

18 

107 

81 

36 

5975.69 

1906 

5241 

2 

22 

7 

17 

146 

804 

37 

5524.00 

1296 

3222 

2 

12 

6 

16 

98 

681 

38 

6222.69 

268 

1050 

3 

12 

6 

16 

60 

1058 

39 

5117.60 

1918 

9662 

4 

12 

7 

23 

120 

1222 

40 

6044.67 

541 

2865 

4 

17 

8 

19 

149 

1409 

* 41 

9126.55 

299 

1903 

5 

20 

7 

35 

186 

1797 

42 

4759.74 

1052 

4280 

3 

12 

4 

20 

77 

653 

43 

6196.51 

1524 

1296 

1 

12 

5 

14 

61 

317 

44 

4602.18 

82 

132 

1 

12 

6 

13 

67 

1800 

45 

4732.80 

64 

309 

3 

9 

4 

15 

51 

777 

* 46 

10335.26 

455 

1316 

2 

16 

7 

2 

91 

1300 

* 47 

11036.11 

1 

9 

9 

0 

0 

0 

1 

0.67 

* 48 

12070.11 

25 

179 

6 

12 

8 

4 

3 

300 

* 49 

8684.86 

1 

11 

11 

0 

0 

0 

1 

0.76 

50 

5926.72 

1476 

5724 

3 

12 

6 

18 

115 

1010 

* tolerance  is  increased  to  5% 


Table  5.  B&B  results  with  0.5%  tolerance  and  discount 


rate  of  0.02 


Problem 

Optimum 

Solution 

Totl.  no 
subprob. 

Total  no 
piv/prob 

Ave  no 
piv/subpr 

max  confl 
set 

Ave.  no 
children 

No  of 
preced 
added/pr 

Max  list 
size 

CPU 

time 

(sec) 

I 

27215.38 

27 

85 

2 

6 

5 

3 

5 

266 

2 

16577.18 

3 

8 

2 

2 

3 

2 

2 

0.42 

3 

21135.17 

5 

4 

4 

4 

5 

1 

2 

0.51 

4 

22493.23 

1 

8 

8 

0 

0 

0 

1 

0.62 

5 

1204.88 

1 

11 

11 

0 

0 

0 

1 

0.76 

6 

3955.70 

13 

100 

7 

12 

13 

2 

13 

6.40 

7 

21053.30 

15 

34 

1 

6 

5 

2 

7 

1.54 

8 

7021.65 

12 

48 

2 

3 

2 

3 

7 

1.65 

9 

6644.25 

957 

5075 

4 

18 

8 

5 

52 

240 

10 

9137.77 

1 

3 

3 

0 

0 

0 

1 

0.21 

81 

Table  5 — continued 


Problem 

Optimum 

Solution 

Totl.  no 
subprob . 

Total  no 
piv/prob 

Ave  no 
piv/subpr 

max  confl 
set 

Ave.  no 
children 

No  of 
preced 
added/pr 

Max  list 
size 

CPU 

time 

(sec) 

11 

335.38 

3 

4 

1 

2 

3 

2 

3 

0.40 

12 

4602.07 

71 

287 

3 

15 

7 

10 

2 

24 

13 

2418.36 

1963 

14958 

6 

19 

5 

5 

48 

983 

14 

5460.82 

736 

4635 

5 

20 

8 

3 

62 

1152 

15 

1322.51 

159 

1334 

7 

12 

6 

3 

29 

331 

16 

2496.04 

405 

3003 

6 

8 

6 

5 

32 

193 

17 

1750.12 

1581 

6248 

3 

12 

7 

8 

69 

615 

18 

2647.25 

109 

771 

5 

6 

3 

4 

25 

51 

19 

2934.07 

106 

356 

2 

12 

6 

17 

73 

331 

20 

1345.65 

35 

158 

3 

8 

5 

7 

26 

693 

21 

3066.95 

1763 

8418 

3 

9 

2 

17 

26 

749 

22 

3344.78 

67 

345 

4 

6 

6 

2 

27 

24 

23 

4366.41 

39 

45 

1 

7 

3 

13 

29 

344 

24 

3526.14 

1077 

2655 

2 

17 

4 

21 

65 

459 

25 

4823.85 

737 

2327 

2 

15 

9 

14 

94 

560 

26 

4348.50 

1755 

4353 

2 

12 

9 

6 

51 

460 

27 

3473.57 

1791 

6582 

3 

20 

4 

18 

101 

615 

28 

4875.39 

1895 

4123 

1 

12 

6 

9 

75 

438 

29 

3506.47 

575 

2607 

3 

8 

3 

5 

34 

203 

30 

2916.13 

115 

421 

3 

7 

4 

5 

21 

341 

31 

5510.54 

128 

576 

4 

16 

8 

17 

109 

1153 

32 

5684.76 

75 

68 

1 

8 

6 

13 

64 

350 

33 

7785.71 

1505 

9639 

5 

12 

6 

13 

89 

1475 

34 

7653.12 

241 

624 

2 

12 

6 

13 

75 

748 

35 

3812.04 

125 

512 

3 

12 

7 

18 

107 

84 

36 

5009.00 

174 

387 

1 

12 

7 

25 

149 

474 

37 

4550.18 

2009 

4179 

1 

12 

6 

16 

98 

691 

38 

5397.56 

1036 

3797 

3 

19 

7 

5 

16 

480 

39 

4422.51 

113 

669 

5 

12 

7 

16 

83 

1218 

40 

5356.67 

511 

2756 

4 

17 

8 

19 

149 

1403 

*41 

8251.99 

13 

46 

3 

12 

13 

1 

2 

43 

42 

4223.62 

958 

3834 

3 

12 

4 

20 

77 

651 

Table  5 — continued 
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Problem 

Optimum 

Solution 

Totl.  no 
subprob. 

Total  no 
piv/prob 

Ave  no 
piv/subpr 

max  confl 
set 

Ave.  no 
children 

No  of 
preced 
added/pr 

Max  list 
size 

CPU 

time 

(sec) 

43 

5420.88 

1358 

1157 

1 

12 

5 

14 

61 

323 

44 

3809.96 

82 

132 

1 

12 

6 

13 

67 

300 

45 

4243.34 

66 

308 

3 

9 

4 

16 

52 

762 

* 46 

8218.00 

95 

281 

2 

12 

6 

16 

82 

1800 

*47 

8213.21 

1 

9 

9 

0 

0 

0 

1 

0.69 

* 48 

9837.96 

143 

1069 

6 

12 

7 

22 

100 

1800 

* 49 

6102.30 

1 

11 

11 

0 

0 

0 

1 

0.76 

50 

5090.03 

1314 

5124 

3 

12 

6 

18 

115 

994 

* tolerance  is  increased  to  5% 


Table  6.  B&B  results  with  0.02  discount  rate 


Problem 
no  of 
adivit. 

Range  of 

Subprob 

solved 

Ave.  no. 
subprob 
per  prob 

Ave.  no. 
pivot/subp 

Ave.  no. 
pivot/prob 

Max. 

conflict 

set 

Ave.  no. 
children 

No.  of 
precedence 
tel.  added 

Ave.  list 
size 

Max.  list 
size 

Ave. 

CPU 

lime 

<22 

[1,9571 

127 

2 

657 

[0,18] 

4 

3 

10 

[1.52] 

30s 

22,23 

[1,2000] 

715 

4 

2729 

[0,20] 

5 

10 

46 

[1,101] 

365s 

27 

[75,2000] 

646 

3 

2529 

[8,19] 

6 

14 

86 

[16,149] 

921s 

35 

[159,736] 

447 

6 

2984 

[12,30] 

7 

3 

45 

[29,62] 

726s 

51 

[13,143] 

251 

5 

465 

12 

7 

[1,22] 

61 

[2,100] 

[43, 

1800] 

Table  7 . B&B  results  with  0.01  discount  rate 


Problem 
no  of 
activit. 

Range  of 

Subprob 

solved 

Ave.  no. 
subprob 
per  prob 

Ave.  no. 
pivot/subp 

Ave.  no. 
pivot/prob 

Max. 

conflict 

set 

Ave.  no. 
children 

No.  of 
precedence 
rel.  added 

Ave.  list 
size 

Max.  list 
size 

Ave. 

CPU 

time 

<22 

[1,1055] 

151 

2 

741 

[0,24] 

4 

3 

14 

[1.63] 

67s 

22,23 

[1,2000] 

776 

4 

2901 

[0,20] 

5 

9 

48 

[1.99] 

364s 

27 

[67,1918] 

861 

3 

3558 

[8,22] 

7 

15 

103 

[57,149] 

873s 

35 

[181,442] 

311 

6 

2123 

[19,20] 

7 

4 

33 

[31,35] 

531s 

51 

[25,455] 

259 

4 

1132 

[12,20] 

7 

[2,35] 

93 

[3,186] 

[310, 

1800] 
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3.3.6  Conclusion 

In  this  section,  we  presented  a depth-first  branch-and- 
bound  algorithm  for  the  resource  constraint  project  scheduling 
problem  with  discounted  cash  flows.  To  our  knowledge  this  is 
the  first  attempt  to  solve  RCPSPDC  in  project  scheduling 
literature  optimally.  The  algorithm  finds  an  optimal  solution 
in  a finite  number  of  steps.  We  tested  the  algorithm  on  50 
test  problems  with  number  of  activities  ranging  from  7 to  51. 
The  following  encouraging  results  were  obtained:  For  problems 
with  less  than  51  activities,  the  algorithm  determined  the 
optimal  solution  with  reasonable  CPU  times.  For  51  activity 
problems,  algorithm  terminated  with  the  best  lower  bounds 
found  within  30  minute  CPU  time  limit.  The  branching  rule 
that  is  based  on  minimal  delaying  alternatives  seems  to  be 
effective  in  bounding  the  size  of  the  tree.  Among  all  the 
problems  solved,  the  maximum  number  of  delaying  alternatives 
is  24  and  average  number  of  branches  generated  from  a node  is 
7 . Branching  from  the  node  with  the  highest  upper  bound 
proved  to  be  effective  in  fathoming  upon  backtracking.  When 
all  the  test  problems  are  considered,  the  maximum  number  of 
unfathomed  nodes  at  any  time  of  the  procedure  is  186  (maximum 
list  size) . 

We  complete  this  section  by  concluding  that  the  branch 
and  bound  procedure  is  a practical  and  efficient  procedure  for 
RCPSPDC  that  produce  optimal  solutions  in  a reasonable 
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computational  effort.  Furthermore,  we  would  like  to  mention 
that  this  is  the  first  exact  solution  technique  for  RCPSPDC. 

3 . 4 Concluding  Remarks 

In  this  chapter,  a heuristic,  Tabu  Search,  and  an  exact 
solution,  Branch  and  Bound,  procedure  for  the  RCPSPDC  were 
developed.  The  computational  results  of  the  Tabu  Search 
procedure  and  Branch  and  Bound  procedure  were  reported.  It 
was  concluded  that  Tabu  Search  and  Branch  and  Bound  procedures 
are  viable  approaches  for  this  problem.  Tabu  Search  produces 
good  near  optimal  solutions  and  Branch  and  Bound  procedure 
yields  optimum  solutions  with  reasonable  computational  effort. 
In  the  next  chapter  we  will  introduce  another  project 
scheduling  problem  which  has  not  been  researched  in  the 
literature.  It  is  named  the  Resource  Constrained  Time/Cost 
Tradeoff  Project  Scheduling  Problem  with  Discounted  Cash 


Flows. 


CHAPTER  4 


THE  RESOURCE  CONSTRAINED  TIME/COST  TRADEOFF  PROJECT  SCHEDULING 
PROBLEM  WITH  DISCOUNTED  CASH  FLOWS 

4.1  Introduction 

In  the  previous  chapter,  we  have  developed  solution 
techniques,  exact  and  heuristic,  for  the  projects  with  limited 
resources  and  precedence  relations.  In  this  chapter,  we  will 
look  into  an  untouched  yet  very  practical  project  scheduling 
problem  named  as  The  Resource  Constrained  Time/Cost  Tradeoff 
Project  Scheduling  Problem  with  Discounted  Cash  Flows, 
RCTCTPSPDC . In  many  real  life  project  scheduling  problems,  the 
project  manager  who  tries  to  determine  the  schedules  that 
maximize  the  net  present  value  of  cash  flows,  has  the  choice 
to  reduce  certain  activity  durations  by  allocating  more 
resources  to  these  activities.  This  reduction  in  the  activity 
durations  (crashing)  may  increase  project  cash  outflow,  but  if 
correct  activities  are  crashed,  gain  in  net  present  value  may 
also  result.  Thus,  the  project  manager  faces  a scheduling 
problem  which  involves  resource  limitations,  variable  activity 
durations  and  precedence  relations.  To  our  knowledge,  no 
research  has  been  done  that  suggests  solution  procedures  for 
this  problem.  In  the  following  section  the  problem 
formulation  will  be  introduced.  In  the  succeeding  sections 
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heuristic  and  exact  solution  procedures  and  their 
implementations  will  be  studied. 

4.2  The  Problem  Formulation 

The  resource  constrained  time/cost  tradeoff  project 
scheduling  problem  with  discounted  cash  flows  involves 
scheduling  and  determining  the  duration  of  project  activities 
with  cash  inflows  and  outflows  in  such  a way  that  the  resource 
and  precedence  constraints  are  met  and  the  net  present  value 
of  the  cash  flows  is  maximized. 


N -*i  Ki  T 

(PR)  Max  £ £ ( Yl  NFik  Yiky>  xite~at-  E ( t-due)  XNtPe-at 

i-1  t-ei  Jc=l  t -due 


subject  to 


E yit=1  Vi 

Jc=l 


E = 1 


N Ki  t+di*-l 

EE(  E ri*vkxiq)  Yik±R. 

i= 1 k=l  q=t 


Vt,  Vav 
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Y±k  e{0, 1}  Vi,  k 


Xit  e(0,  l)  Vi,  t 


i. 


E 

t = (=> 


tXmC  + y£  d nkYnkZ  £ tXBt 


k=l 


t=e„ 


( m,n ) e H 


where 

T : An  upper  bound  on  the  completion  time  of  the  project, 
due  : Due  date  on  the  completion  time  of  the  project. 

P : Penalty  per  period  of  delaying  the  project  beyond  its 

due  date. 

N : Last  activity  in  the  project. 

qt.  : Cash  received  at  the  completion  time  of  activity  i. 

K1  : Number  of  possible  durations  for  activity  i. 
d^  : kth*  possible  duration  of  activity  i. 

riavk  : Numker  of  units  of  resource  av  required  by  activity  i 
that  has  djk  duration. 

Rav  : Availability  level  of  resource  av,  av=l,...,V. 

NF^  : Net  cash  flow  at  the  completion  time  of 
activity  i that  has  dt.k  duration. 
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NF 


s-vt  + r,  +0^ 


Y; 


k 


Cash  flow  for  activity  i in  period  j,  j=l,  . . . ,djk. 
Total  crashing  cost  if  d..k  duration  is  chosen  for 
activity  i. 

1,  If  djk  duration  is  chosen  for  activity  i; 
,0,  otherwise. 


ri,  If  activity  i is  completed  at  time  t; 

it  • 


. 0,  otherwise. 
a : Discount  rate. 

ei  : Earliest  completion  time  for  activity  i. 
lj  : Latest  completion  time  for  activity  i. 

As  explained  in  section  3.1,  in  order  to  reduce  the 
number  of  integer  variables,  CPM  method  is  used  to  determine 
the  earliest  and  the  latest  activity  completion  times. 

The  first  set  of  constraints  insures  that  only  one 
duration  alternative  is  chosen  among  the  set  of  Ki  possible 
durations  for  each  activity  of  the  project.  The  second  set  of 
constraints  requires  that  each  activity  is  completed  exactly 
once.  The  third  set  of  constraints  represents  the  resource 


89 


constraints.  At  any  time  period  t,  the  total  resource 
consumption  of  the  activities  which  are  scheduled  to  be  active 
at  time  t,  can  not  exceed  the  resource  availability  at  time  t 
of  the  resource  type  av.  Note  that  the  resource  requirement 
of  an  activity  varies  with  the  duration  chosen  for  that 
activity.  The  forth  set  of  constraints  represents  the 
precedence  relations.  If  activity  m immediately  precedes 
activity  n that  has  a duration  dnk,  then  the  completion  time 
of  activity  m plus  dnk  must  be  less  than  or  equal  to  the 
completion  time  of  activity  n.  The  fifth  and  sixth  sets  of 
constraints  are  the  integrality  requirements  on  the  variables 
Y.k  and  Xit  respectively. 

The  first  term  in  the  objective  function  is  the  sum  of 
the  values  of  cash  flows  (including  crashing  cost)  discounted 
to  time  zero.  The  second  term  is  the  discounted  value  of 
delay  penalties  that  occurs  if  the  project  duration  exceeds 
its  due  date. 

If  the  variable  duration  assumption  is  relaxed  in  (PR) , 
then  the  problem  is  a RCPSPDC.  Thus,  the  procedures  that 
determine  a due  date  and  the  latest  completion  time  T of  a 
project  in  section  3.2.2  can  be  used  to  determine  a due  date 
and  T for  (PR)  problem. 

In  the  model,  the  following  assumptions  are  made: 

1.  Each  activity  is  assumed  to 

- have  Ki  possible  integer  durations; 
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- require  known  level  of  resources  of  each  type  which 

increases  as  the  duration  of  an  activity  shortens; 

- utilize  constant  level  of  each  resource  type  throughout 

its  duration; 

- have  a linear  crashing  cost  curve; 

- have  no  preemption. 

2.  The  level  of  availability  of  each  resource  type  is  assumed 
to  be  constant  throughout  the  schedule. 

4.3  A Heuristic  Procedure  for  RCTCTPSPDC 

In  this  section,  we  present  a heuristic  procedure  for  the 
resource  constrained  time/cost  tradeoff  project  scheduling 
problems  with  discounted  cash  flows.  The  logic  of  the 
procedure  is  similar  to  that  of  the  one  presented  in  section 
3.2.2.  The  procedure  performs  a single  forward  pass  through 
the  network  starting  with  the  first  activity  which  is  usually 
a dummy  activity.  At  the  completion  time  of  an  activity,  all 
the  activities  which  can  be  scheduled  according  to  their 
precedence  relations  are  determined.  Let  L be  the  set  which 
contains  the  schedulable  activities.  The  best  duration 
alternative  is  determined  for  each  activity  in  L by  computing 
an  index  called  the  Net  Marginal  Gain  (NMG) . It  gives  the 
change  in  the  net  present  value  in  one  time  unit  per  unit 
increase  in  the  resource  usage  when  different  duration 
alternatives  are  considered  for  an  activity.  The  durations 
with  the  maximum  NMG  values  are  assigned  to  the  activities  in 
L.  Let  M be  the  sorted  set  of  activities  in  L according  to 
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their  NMG  values.  The  procedure  begins  scheduling  the 
activities  in  M until  a resource  conflict  occurs.  At  this 
instant,  unscheduled  activities  in  M are  delayed  and  procedure 
moves  to  the  next  activity  completion  time.  The  procedure 
continues  in  this  fashion  until  the  last  activity  which  is 
also  a dummy  activity  is  scheduled. 


4.3.1.  A Formal  Statement  of  the  Procedure: 

1.  Start  with  scheduling  the  initial  project  activity  usually 
a dummy  activity  at  t=0. 

2.  At  the  completion  time,  t,  of  an  activity,  only  considering 
the  precedence  relations,  determine  the  set  L which  contains 
all  the  schedulable  activities  that  are  currently  at  their 
normal  durations. 

3.  If  L=0,  THEN  move  to  the  next  completion  time. 

If  the  last  scheduled  activity  is  the  last 
activity  THEN  return  the  schedule  and  the  net  present 
value  and  STOP. 

ELSE  go  to  step  2 . 

ELSE  go  to  step  4 . 

4.  For  each  activity  i e L,  determine. 


-o  t, 


- 


* _ g.  ( e~g Cj  -e-gtM 


to  ties 
j- 1 


k = l , . . , K± . 


where 


92 

H;k  : the  net  marginal  gain  of  activity  i which  has 
duration. 

tjk  : the  completion  time  of  an  activity  i which  has  d^ 
duration. 

tj  : the  completion  time  of  activity  i that  has  normal 
duration. 

5.  Determine  the  best  duration  for  each  activity  in  L. 


Hk'  = max  {H±  k : k = 1 , . . . } 


Let  d{k*  be  the  duration  chosen  for  activity  i; 

6.  Let  M be  the  list  of  sorted  activities  in  decreasing  order 
of  their  H^*  values. 

7.  Schedule  the  activities  at  time  t starting  with  the  first 
element  of  M until  a resource  conflict  occurs. 

7.1.  If  no  resource  conflict  occurs  THEN  schedule  all  the 
activities  in  M and  go  to  step  8, 

ELSE  if  a resource  conflict  occurred  during  the 
scheduling  of  activity  j , then  activity  j and  all  the 
succeeding  elements  of  M are  delayed,  go  to  step  8. 

8 . Go  to  step  2 . 

As  it  is  clear  from  step  7 of  the  procedure,  the  schedule 
time  of  each  activity  depends  solely  on  its  NMG  value.  The 
bigger  the  NMG  value  is,  the  better  the  chance  of  an  activity 
to  be  scheduled  earlier.  The  NMG  value  of  an  activity 
furthermore  depends  on  resource  utilization  and  crashing  cost 
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of  that  activity  when  d^*  duration  is  chosen.  Thus,  high 
resource  requirements  and/or  costly  crashings  and/or  negative 
cash  flows  may  result  in  delays  in  activity  schedule  times. 

Note  that  this  procedure  will  always  produce  a feasible 
solution  to  the  RCTCTPSPDC. 

4.3.2.  Computational  Experience 

The  above  heuristic  procedure  was  coded  in  Fortran  and 
tested  on  the  50  test  problems  used  in  sections  3.2.5  and 
3.3.5.  The  computational  experiments  were  done  on  an  IBM 
3090-600J  computer.  The  following  data  were  added  to  the 
original  data  set  prepared  for  the  50  problems.  For  each 
activity  three  durations  are  considered;  minimum,  average  and 
normal  durations.  The  minimum  and  average  durations  are 
obtained  by  reducing  the  normal  activity  durations  by  60%  and 
40%  respectively.  If  an  activity  has  a normal  duration  of  1 
time  unit  then  minimum  and  average  durations  are  set  to  the 
normal  duration.  Similarly,  the  resource  requirements  for  the 
activities  at  their  average  or  minimum  durations  are  increased 
to  50%  and  80%  of  the  normal  resource  requirements 
respectively.  Furthermore,  if  the  normal  resource  requirement 
of  an  activity  is  zero,  then  the  average  and  the  minimum 
resource  requirements  are  set  to  1 and  2 respectively.  The 
crashing  costs  per  unit  time  are  determined  as  follows:  Cash 
received  at  the  completion  time  of  the  activities  are 
generated  randomly  and  uniformly  between  [-500,1000]  in  the 
original  data  set,  (see  section  3.2.5).  Thus,  at  the  earliest 
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completion  time  of  any  activity,  which  is  t=l,  the  maximum 
change  in  the  net  present  value  is  (1000-1000*exp (-0 . 01*1) ) = 
10  with  a=0 . 01  and  ( 1000-1000*exp (-0 . 02*1)  = 20  with  a=0.02. 
Therefore,  to  make  a tradeoff  between  cash  received  and  the 
cost  of  reduction  in  the  activity  time,  the  crashing  cost  for 
each  activity  per  unit  time  was  generated  randomly  between 
[1,10].  Note  that,  it  was  assumed  that  the  cash  flows,  f^s, 
are  negligible. 

The  computational  results  are  given  in  table  8.  The  net 
present  values  are  given  in  column  2.  The  first  row  in  each 
cell  corresponds  to  the  results  obtained  by  using  the  discount 
rate  of  0.01  and  the  second  row  corresponds  to  the  results 
obtained  by  using  the  discount  rate  of  0.02.  Column  3 gives 
the  total  number  of  activities  crashed  to  their  average 
durations  and  column  4 gives  the  total  number  of  activities 
crashed  to  their  minimum  durations.  Last  column  gives  the  CPU 
time  that  each  network  consumed  in  seconds. 

In  all  the  problems,  the  total  number  of  activities 
crashed  to  their  average  durations  is  greater  than  or  equal 
to  the  number  of  activities  crashed  to  their  minimum 
durations.  Since  the  maximum  crashing  cost  per  time  unit  is 
set  to  10,  as  one  would  suspect,  the  total  number  of 
activities  crashed  from  their  normal  durations  with  0.02 
discount  rate  for  each  problem  is  greater  than  or  equal  to  the 
total  number  of  activities  crashed  from  their  normal  durations 


with  0.01  discount  rate. 


95 


Table  8 . Computational  results  obtained  by  the  heuristic 
procedure 


Problem 

Heuris. 

Soln. 

No.  of 
act.  AVE 
crashed. 

No.  of 
act.  min 
crashed 

CPU 

Times 

(sec) 

1 

30150.17 

0 

0 

0.11 

27101.17 

0 

0 

0.11 

2 

17213 . 01 

2 

0 

0.11 

16705.23 

2 

0 

0.11 

3 

21510.37 

6 

0 

0.11 

18737.10 

7 

0 

0.12 

4 

22735.08 

0 

0 

0.12 

21486.04 

0 

0 

0.14 

5 

1049.94 

0 

0 

0.14 

1004.00 

0 

0 

0.14 

6 

3965.77 

0 

0 

0.15 

3664.04 

0 

0 

0.14 

7 

22398.25 

4 

0 

0.15 

21849.72 

5 

0 

0.11 

8 

7522 . 07 

2 

0 

0.11 

6844.85 

2 

0 

0.11 

9 

6797.75 

6 

0 

0.11 

5792.60 

7 

0 

0.14 

10 

10545.50 

2 

2 

0.11 

9984.01 

2 

2 

0.11 

11 

356.06 

2 

0 

0.12 

306.53 

2 

0 

0.11 

12 

4766.78 

2 

0 

0.12 

4440.55 

5 

0 

0.15 

13 

2220.66 

5 

1 

0.17 

2001.44 

6 

1 

0.16 

14 

6654.85 

8 

0 

0.23 

5560.71 

11 

0 

0.24 

15 

1679.53 

7 

0 

0.26 

1042.68 

12 

0 

0.25 

16 

2291.60 

2 

0 

0.16 

1878.95 

4 

1 

0.16 

Table  8 — continued 


Problem 

Heuris. 

No.  of 

No.  of 

CPU 

Soln. 

act.  AVE 

act.  min 

Times 

crashed. 

crashed 

(sec) 

17 

2013 . 98 

2 

0 

0.16 

1590.44 

1 

0 

0.16 

18 

2783.08 

7 

0 

0.17 

2459.27 

7 

0 

0.16 

19 

3573.49 

5 

1 

0.20 

2968.62 

6 

2 

0.16 

20 

1654.58 

3 

1 

0.16 

1345.58 

4 

1 

0.15 

21 

3848 .36 

5 

2 

0.16 

3209.34 

6 

3 

0.16 

22 

3579.32 

5 

0 

0.17 

2961.67 

6 

1 

0.16 

23 

5197.10 

3 

0 

0.16 

4210.90 

8 

0 

0.16 

24 

4612.62 

8 

0 

0.16 

3655.90 

10 

0 

0.16 

25 

5656.21 

7 

0 

0.16 

4724.80 

11 

0 

0.15 

26 

4773.47 

5 

1 

0.16 

4196.13 

6 

1 

0.16 

27 

4099.27 

5 

1 

0.16 

3462.69 

6 

1 

0.16 

28 

5181.33 

6 

0 

0.16 

44517.42 

7 

0 

0.15 

29 

3942.25 

5 

1 

0.16 

3581.20 

6 

1 

0.16 

30 

3222.62 

7 

0 

0.16 

2668.88 

9 

0 

0.16 

31 

6623.73 

7 

1 

0.19 

5434.11 

12 

1 

0.18 

32 

6749.51 

9 

0 

0.19 

5769.53 

14 

0 

0.18 

33 

8584.99 

13 

0 

0.19 

7334.73 

14 

1 

0.18 

Table  8 — continued 


Problem 

Heuris. 

No.  of 

No.  of 

CPU 

Soln. 

act.  AVE 

act.  min 

Times 

crashed. 

crashed 

(sec) 

34 

8433.51 

8 

0 

0.19 

7170.75 

12 

0 

0.18 

35 

4594.32 

10 

0 

0.18 

3524.28 

14 

0 

0.18 

36 

5994.22 

7 

0 

0.18 

4931.56 

9 

0 

0.18 

37 

5127.34 

5 

1 

0.18 

4263.14 

9 

0 

0.18 

38 

6012.03 

7 

0 

0.18 

5050.33 

10 

0 

0.18 

39 

4925.65 

5 

2 

0.19 

3803.07 

9 

2 

0.18 

40 

5731.48 

8 

0 

0.18 

4931.95 

11 

0 

0.18 

41 

8801.58 

13 

1 

0.49 

6299.41 

21 

0 

0.46 

42 

4633.04 

5 

1 

0.16 

4214.35 

7 

1 

0.16 

43 

6008.01 

9 

0 

0.16 

5171.61 

11 

0 

0.15 

44 

4380.30 

3 

0 

0.16 

3313.08 

6 

0 

0.16 

45 

4544.57 

4 

2 

0.16 

4239.79 

8 

2 

0.16 

46 

8845.07 

9 

1 

0.49 

6507.77 

12 

2 

0.45 

47 

10533.08 

17 

1 

0.49 

7824.08 

28 

0 

0.44 

48 

11530.12 

11 

0 

0.47 

8474.60 

18 

0 

0.43 

49 

8345.92 

12 

0 

0.45 

6233.56 

23 

0 

0.44 

50 

5639.65 

9 

0 

0.19 

4615.06 

11 

0 

0.18 
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As  it  is  clear  from  step  5 of  the  algorithm,  the  best 
duration  chosen  for  each  activity  depends  on  the  way  the 
resource  consumption  and  the  crashing  cost  values  are  chosen. 
Therefore,  to  show  the  effect  of  the  values  of  these 
parameters  on  the  solution  quality,  the  following  changes  were 
made:  the  resource  consumption  of  an  activity  with  average  and 
minimum  durations  were  increased  to  40%  and  50%  of  the  normal 
resource  requirements.  Furthermore,  crashing  costs  were 
generated  randomly  and  uniformly  between  [1,5]  and  between 
[1,20].  Each  of  these  changes  were  tested  on  the  arbitrarily 
chosen  10  test  problems  among  the  50  test  problems.  The 
following  results  were  obtained: 

Table  9. Results  with  40%  and  50%  increase  in  resource 


consumption. (a=0.02) 


Problem 

Heuristic 

Solution 

Ave . Crsh 
activities 

Min. crsh 
activit . 

CPU  time 
(sec. ) 

2 

16705.23 

2 

0 

0.13 

3 

18935.01 

5 

2 

0.13 

4 

21936.90 

0 

0 

0.15 

5 

1013.80 

0 

0 

0.13 

50 

4767.89 

7 

4 

0.19 

7 

21849.72 

5 

0 

0.12 

19 

3191.13 

2 

6 

0.18 

20 

1345.58 

4 

1 

0.18 

22 

3131.57 

5 

2 

0.18 

47 

8657.40 

26 

2 

0.44 
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Table  10. Results  with  40%  and  50%  increase  in  resource 


consumption  and  crashing  cost  is  between  [1,5] 
(g=0. 02) . 


Problem 

Heuristic 

Solution 

Ave . Crsh 
activities 

Min. crsh 
activit . 

CPU  time 
(sec. ) 

10 

10015.58 

2 

2 

0.12 

11 

338.15 

3 

0 

0.12 

12 

4489.72 

5 

0 

0.16 

13 

2000.80 

7 

1 

0.18 

14 

6087.11 

13 

5 

0.24 

15 

1152.11 

12 

2 

0.26 

35 

3639.19 

11 

3 

0.19 

36 

5542.66 

5 

7 

0.19 

37 

4447.01 

9 

2 

0.20 

38 

5316.48 

11 

2 

0.20 

Table  11. Results  with  40%  and  50%  increase  in  resource 


consumption  and  crashing  cost  is  between  [1,20] 
(g=0. 01) . 


Problem 

Heuristic 

Solution 

Ave . Crsh 
activities 

Min. crsh 
activit. 

CPU  time 
(sec. ) 

40 

5668.34 

5 

3 

0.20 

41 

9086.56 

4 

3 

0.48 

42 

4613.70 

3 

1 

0.17 

43 

6058.28 

5 

2 

0.17 

44 

4378.18 

3 

0 

0.17 

45 

4465.23 

2 

1 

0.17 

46 

9565.21 

6 

1 

0.47 

47 

8657.40 

26 

2 

0.47 

48 

11539.50 

7 

0 

0.47 

49 

8555.92 

4 

0 

0.46 
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Table  12. Results  with  50%  and  80%  increase  in  resource 


consumption  and  crashing  cost  is  between  [1,20] 
(a=0. 01) . 


Problem 

Heuristic 

Solution 

Ave.  Crsh 
activities 

Min. crsh 
activit . 

CPU  time 
(sec. ) 

40 

5659.07 

5 

0 

0.20 

41 

8741.08 

7 

0 

0.48 

42 

4613.70 

3 

1 

0.17 

43 

5982 . 62 

7 

0 

0.17 

44 

4377.41 

3 

0 

0.17 

45 

4464.67 

3 

0 

0.17 

46 

9317.50 

6 

1 

0.47 

47 

10533.73 

10 

0 

0.47 

48 

11426.71 

7 

0 

0.47 

49 

8251.24 

4 

0 

0.46 

In  general,  reduction  in  the  resource  consumption  alone 
increased  the  number  of  activities  crashed  to  their  minimum 
durations  and  resulted  in  an  increase  in  the  optimum  solution 
values.  Reduction  in  the  crashing  costs  along  with  the 
reduction  in  the  resource  consumptions  resulted  in  an  increase 
in  the  optimal  solution  values  and  the  number  of  activities 
crashed  to  their  minimum  and  average  durations.  Finally,  an 
increase  in  the  crashing  costs  resulted  in  a decrease  in  the 
optimal  solution  values.  The  following  table  gives  the 
percentage  changes  in  the  optimal  solution  values  and  in  the 
number  of  activities  crashed. 
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Table  13. 


A 

B 

C 

D 

Average 

18% 

8% 

28% 

36% 

Duration 

increase 

increase 

decrease 

decrease 

Minimum 

325% 

87% 

66% 

66% 

Duration 

increase 

increase 

increase 

decrease 

Heuristic 

70% 

90% 

50% 

80% 

Solution 

increase 

increase 

decrease 

decrease 

A:  a=0.02,  40%  and  50%  increase  in  resource  consumptions, 
crashing  costs  are  generated  between  [1,10]. 

B:  a=0.02,  40%and  50%  increase  in  resource  consumptions, 
crashing  costs  are  generated  between  [1,5]. 

C:  a=0.01,  40%  and  50%  increase  in  resource  consumptions, 
crashing  costs  are  generated  between  [1,20]. 

D:  a=0.01,  50%  and  80%  increase  in  resource  consumptions, 
crashing  costs  are  generated  between  [1,20]. 

4.3.3.  Conclusions 

In  this  section,  a single  pass  heuristic  procedure  is 
developed  for  the  resource  constrained  time/cost  tradeoff 
scheduling  problems  with  discounted  cash  flows.  The 
performance  of  the  procedure  was  tested  on  the  50  test 
problems.  Computational  testing  indicated  that  the  procedure 
produces  feasible  solutions  with  modest  computational  effort. 
Furthermore,  the  effects  of  assigning  different  values  to 
parameters  are  analyzed  on  arbitrarily  chosen  sets  of 
problems.  In  general,  it  was  concluded  that  reducing  the 
resource  requirement  of  the  activities  and  reducing  the 
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crashing  cost  per  unit  time  resulted  in  an  increase  in  the 
objective  function  values. 

4.4  A Subgradient  Optimization  Approach  to  RCTCTPSPDC 
4.4.1  Introduction 

In  this  section,  we  will  develop  a Subgradient 
Optimization  algorithm  to  solve  RCTCTPSPDC  optimaly. 
Subgradient  Optimization  technique  is  commonly  used  for 
solving  the  Lagrangian  dual  problem  of  a nonlinear  programming 
problem  [3].  The  following  definitions  about  a Lagrangian 
duality  will  be  needed  to  facilitate  a background  in  this 
technique. 

The  Lagrangian  Dual  Problem: 

Let  the  primal  problem  (PP)  be: 

(PP)  Maximize  f(x) 

Subject  to 

g (x)  < b 
x e X 

Then,  the  Lagrangian  Dual  problem  of  problem  (PP)  is  (LD) : 
(LD)  Minimize  0(u,v) 

Subject  to  u > 0 

where  0(u,v)  = max  {f(x)+  u*  (b-g(x))  :x  e X). 

The  Lagrangian  Dual  problem  is  obtained  by  incorporating 
the  constraint  set  g(x)  in  the  objective  function  with 
corresponding  multipliers,  u,  called  Lagrangian  multipliers. 
A problem  might  have  several  different  relaxations  and 
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accordingly  there  might  be  several  different  approaches  for 
solving  these  relaxations.  The  subgradient  optimization 
technigue  uses  the  idea  that  moving  a small  step  in  a 
subgradient  direction  gets  us  closer  to  an  optimal  dual 
solution.  Under  certain  convexity  assumptions,  the  optimal 
dual  solution  obtained  is  an  optimal  solution  to  maximization 
problem  (PP) . The  following  flow  chart  summarizes  this 
technique  for  a maximization  problem  with  one  set  of  less  than 
or  equal  to  type  constraints. 


Figure  4.  An  algorithmic  framework  of  Subgradient  Optimization 
Technique 
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In  the  procedure,  it  is  presumed  that  lagrangian  relaxing 
the  complicated  constraints  would  make  the  resulting  problem 
easier  to  solve.  Thus,  the  efficiency  of  the  algorithm  solely 
depends  on  the  computational  efficiency  of  the  technique  used 
to  solve  the  subproblem  at  each  iteration.  Before  moving  to 
the  next  iteration,  the  lagrangian  multipliers  are  updated  by 
choosing  a suitable  step  size  tk.  There  are  two  commonly  used 
step  sizes.  Both  are  developed  by  Polyak  [61]. 

Step  Size  1. 

tk  = — — 

k II  9k  II 


Step  Size  2. 


ek  (6  ( u k)  -6*) 

II  <7*  II2 


where 

0 < ek  < 2 and  0*  is  a lower  bound  on  0 that  can  be  obtained 
by  applying  a heuristic  to  the  primal  problem  (PP) . In 
general,  empirical  results  indicate  that  step  size  2 yields 
superior  results  compared  to  the  results  obtained  with  other 
step  sizes  even  though  it  requires  a lower  bound  which  might 
be  difficult  to  find  for  some  problems.  For  details  of  step 
sizes  and  the  applications  of  subgradient  optimization 
technique  refer  to  Bazaara  and  Shetty  [3],  Polyak  [61]  and 
Venkataramanan  et  al.  [80]. 
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4.4.2  The  Application  of  Subgradient  Optimization  Technique  to 

RCTCTPSPDC 

In  this  section,  first  a lagrangian  relaxation,  (PR1)  of 
(PR)  will  be  given.  It  is  obtained  by  relaxing  the  resource 
constraints  and  putting  this  set  of  constraints  in  the 
objective  function  with  lagrangian  multiplier  vector,  X.  It 
will  also  be  shown  that  the  relaxed  problem,  (PR'),  can  be 
written  as  a Time/cost  Tradeoff  Project  Scheduling  Problem 
with  Discounted  Cash  flows  (TCTPSPDC) . It  will  be  further 
shown  that,  following  the  work  of  Erenguc,  Tufekci  and  Zappe 
[19],  the  Generalized  Benders  Decomposition  technique  can  be 
used  to  solve  (PR1). 

Let  (PR')  be  the  problem  obtained  by  lagrangian  relaxing 
the  resource  constraints  of  problem  (PR)  and  by  including  the 
penalty  value  of  the  objective  function  as  a cash  outflow  for 
the  last  activity  of  the  network. 

(PR') 


N -L 

Max  E E NFik  Yik)  Xlte-at  + 

l^l  t=ei  k=l 

T Codes  N Ki  C+di*-l 

E E ( E *i»vxiq)  Yik] 

C= 1 av= 1 i=  1 k=l  q=t 


subject  to 


E Yik  = i Vi 


Xit  e(0,  l)  Vi,  t 
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ii 


E 

t=e. 


*it  = l 


Vi 


E 


E ^ 


t=e„ 


(m,n)  eH 


Y±k  e(o, 1}  Vi,  A: 


Note  that,  in  (PR')  the  lagrangian  relaxed  resource 
constraints  are  written  such  that  the  first  sum  is  from  t=l 
to  T.  However,  recall  that  in  (PR) , we  define  a Xjt  variable 
for  each  integral  time  period  t between  the  earliest  e,.  and 
the  latest  1.  completion  time  of  activity  i.  Thus,  the 
lagrangian  relaxed  constraints  in  the  objective  function  of 
(PR')  can  be  rewritten  such  that  needless  notation  can  be 
discarded.  Therefore,  for  fixed  A,  A'  the  objective  function 
is  equivalent  to 

N Kl  ll 

= Max  E E [ E NFikXite-t  + 

i=  1 k=l  t=ei 

Coties  minfij  , t+dj*- 1) 

E ria/<E  E ] V 

av=l  t=e±  q=  t 


Note  that  the  following  constant  term  is  dropped  from  the 
objective  function: 
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Then, 


T Codes 


E E 


t=l  av=  1 


= Max  £ E E [7/F/e-atXit  + 


i=l  Jc=l  t=ei 
toties  , t+d^^-l} 


E E 


av»l 


<j=C 


W *i  ij 

- ^ EE  E 

i=l  Jc=l  C=ej 

Codes  cxinlli , C+d1k~l)  1*  r * 

E E -( 


av=l 


g=t 


Let  the  constant  term  in  the  bracket  be  NF^.  Note  that  NFitk 
denotes  the  net  future  value  of  all  cash  flows  associated  with 
activity  i if  the  kth‘  duration  is  selected  and  the  activity  is 
completed  at  time  period  t.  Then, 

= Max  E E E 

i=l  k= 1 C=e1 


E ^ = 1 


Jc=l 


Vi 


subject  to 
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E *ic  = 1 Vi 


t=e1 


Xlt  e(0,  l)  Vi,  t 


Y±k  e{0, 1}  Vi,  k 


E ^t  + E <*«**»**  E (m,n)eH 


In  a recent  study  by  Erenguc,  Tufekci  and  Zappe  [19],  the 
time/cost  tradeoff  project  scheduling  problem  with  discounted 
cash  flows  is  formulated  as  a mixed  integer  nonlinear  program 
and  was  solved  by  using  the  Generalized  Benders  Decomposition 
(GBD)  technique  developed  by  Geoffrion  [26] . The  problem 
involves  scheduling  the  activities  with  cash  flows  where 
shorter  activity  durations  can  be  obtained  by  incurring  higher 
costs.  The  objective  of  the  problem  is  to  maximize  the  net 
present  value  of  cash  flows. 

The  following  formulation  proposed  by  the  authors  enabled 
the  GBD  technique  to  be  used  to  solve  TCTPSPDC. 


(PI)  <\ 


G 


'J-I[  (1-e'07")  / ( e“  - 1 ) ] 


(i.j)eB(j)  Jc= 1 


subject  to 
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Ku 

Ti-Tj  £ j =2 , . . , M (i,j)  efl(j) 

Jc-l 


TM  ±R 


£ Yuk  = 1 

Jc*l 


(i,j)  e-A 


1^*  G { 0 , 1 } 


Jc— 1 / . . , 

o 


(i, j) 
i 6 N 


where 


)*V 

B ( j ) 

o 

Qj^' 

A 

XR  P 

N 


Ti 


: The  net  future  value  of  activity  (i,j)  if  d1-k 
duration  is  chosen  for  it. 

: All  activities  leading  into  event  j . 

: Indirect  cost  per  unit  time  paid  at  the  end  of  each 
time  period. 

: The  payment  received  with  the  realization  of  event  j. 
: Set  of  activities. 

: Due  date  on  the  project  completion  time. 

: Set  of  events. 

: Occurrence  time  of  event  j . 

: Nominal  interest  rate. 
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Yii 


k . 


1,  if  the  kth‘  duration  is  selected  for  activity  (i,j)  ; 


,0,  otherwise. 

M : Last  activity  of  a network. 

The  above  formulation  (PI)  is  rewritten  [19]  in  a form 
which  is  more  amenable  to  solution  with  GBD  technique.  The 
details  of  the  GBD  technique  and  the  applicability  of  GBD  to 
(PI)  are  given  in  Appendix  B. 


Claim:  If 

X „ <-  0,  if  =0  for  all  t for  some  av;  otherwise 

tav  ' tav 

Atav  <_  mint  (*ta v/e'at  : (Atav/e"at)  >0  > for  SOIne  aV' 

av=l,...,V,  then  problem  (PR')  can  be 
written  as  (PI)  problem. 

In  what  follows  for  notational  convenience,  we  will 
assume,  without  loss  of  generality,  that  I and  QjS  are  zero 
for  all  j . 

Proof  of  Claim:  Given  the  earliest  and  the  latest  occurrence 
time  of  each  event, ei  and  1 j , problem  (PI)  can  be  written  as 
follows : 

(PI ' ) 

£ £ NFljkYijk)  e-atWjc 

j- 1 C=ej  ( i,j ) 6B(j)  k-1 


subject  to 


Ill 


E tWuc+Y,  duvk  Yuvk  * E tP/vC  V 


t-e,. 


Jc*l 


t=ev 


E ^i/  = l V (i,j)  €A 

*=1 


Yy*  e {0,1}  V(i,j)  6A,  VJc 


where 


Wjt  : 


1,  if  event  j is  realized  at  time  t; 


0,  otherwise. 


Furthermore  to  be  consistent,  we  make  the  following 
change  in  the  notation  of  (PR1):  Let  (i,j)  indicate  an 

activity  leading  into  event  j.  Then  (PR')  can  be  rewritten  as 
(PR") 


1ij  KiJ 


Max  E E E^t%te"^: 

(i,j)  e.A  Jc=l 


U 


subject  to 


E Yijk=1  V(i,j)  EA 

k-1 


E E 


)c=l 


t=ev 


(in,  v)  , ( V,  u)  EH 
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Xijt  £ {0,1}  V(i,j)  , t 
Y1:jk  £ (o,  1 } V (i,  j)  , k 

Problem  (PR")  involves  determining  the  completion  time 
and  the  duration  of  each  activity  such  that  the  net  present 
value  of  cash  flows  that  occur  at  the  completion  time  of  the 
activities  is  maximized.  The  first  set  of  constraints  ensures 
that  only  one  duration  alternative  is  chosen  for  each 
activity.  The  second  set  of  constraints  guarantees  that  each 
activity  completes  exactly  once  and  finally  the  third  set  of 
constraints  represents  the  precedence  relations  between  the 
activities . 

Problem  (PI1)  involves  determining  the  occurrence  time  of 
each  event  and  the  duration  of  each  activity  such  that  the  net 
present  value  of  cash  flows  on  the  events  is  maximized  subject 
to  precedence  constraints  and  the  constraints  which  ensure 
that  each  event  occurs  exactly  once. 

But  by  the  result  of  proposition  3.3.1,  we  know  that  the 
cash  flows  on  the  activities  can  be  equivalently  represented 
by  the  cash  flows  on  the  events.  Thus,  the  objective  function 
and  the  constraints  of  (PR")  is  rewritten  in  the  following 


fashion. 
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Max  £ £ [ £ Y^PNF^Y, 

3- 1 U.v)€F(j)  k=l 


k v k~-*djvk 
jv  e 


JC 


i# 


+ E E^JJt***,*  1 e_BtW/t 

(i,i)€fl(j)  Jc-1 


E ^Ut+E  duvkYuvkz  £ t^vc  V (u,v)  e/l 


C=eu 


ic-1 


C=ev 


E *^*  = 1 V (i,j)  eA 

Jc=l 


y,/  6 {0,1}  V(i,j)6A,Vk 


where 

F ( j ) : The  set  of  activities  emanating  from  event  j. 


B ( j ) : The  set  of  activities  leading  into  event  j. 

and  where 


PNF- 


JVt 


NNF,jt‘ 


' NFjvtk  , if  NFjvtk  > 0; 

. 0,  otherwise. 

NFjjtk  , if  NFijtk  < 0; 


, 0,  otherwise. 

Note  that  in  (PI'),  cash  flows  on  the  nodes,  NFijk,  are 
independent  of  event  times  whereas  in  (PR")  cash  f lows , NFi Jtk, 
vary  according  to  the  time  each  event  occurs.  Recall  that 
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Atav/e  at  term  of  NFj Jtk  depends  on  the  event  occurrence  time.  To 
make  NFjjtk  independent  of  t we  now  define  A,tav  as  follows: 

Atav  <-  0,  if  ltav  =0  for  all  t,  for  some  av;  eq.4.4.2 

*tav  <_  mint^tav/e"at  ^tav/e"at)  > °>  for  S01ne  aV ' 

av=l, . . . ,V. 

That  is,  if  we  assign  the  same  multiplier,  determined  as 
above,  to  all  time  periods  for  each  resource  type  then  it 
becomes  clear  that  the  objective  function  of  (PR")  can  be 
written  as  (PI').  Note,  however  that,  this  approximation  may 
result  in  providing  looser  bounds  on  the  objective  function 
values  of  RCTCTPSPDC. 

Note  that,  the  precedence  relationships  of  (PR")  and 
(PI')  are  of  exactly  the  same  form.  Furthermore,  by  the 
argument  given  in  proposition  3.3.1,  the  precedence 
constraints  of  (PI')  can  be  written  as  the  precedence 
constraints  of  (PI).  Hence  (PR"),  therefore  (PR')  can  be 
written  as  (PI')  therefore  (PI).  Q.E.D. 

The  Subqradient  Algorithm 

Let  problem  (PR)  be  written  in  the  following  simplified  form: 
(S)  Maximize  NPV(x,y) 

subject  to  h(x,y)  < a 
g (x, y)  < b 
x eX,  y eY. 

where  h is  the  set  of  precedence  relations  and  g is  the  set  of 
resource  availability  constraints. 

Then,  the  Lagrangian  dual  problem  of  (S)  is 
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(LS)  Minimize  0(A) 

subject  to  A > 0 

where  0(A)  = maxY  u (NPV(x,y)+AT(b-g(x,y)  ) : h(x,y)<a,  xe  X, 

ye  Y }. 

Step  0.  Select  A0  >0.  Determine  a lower  bound,  w.  Let  k=0. 
Step  l.  Solve  the  following  subproblem  by  using  GBD : 

(SA)  Maxxy  {NPV (x, y) +Ak (b-g (x, y) ) : h(x,y)<a,  xeX, 
yeY). 

Let  xk,yk  constitutes  the  optimal  solution. 

Step  3.  Determine  0(Ak)=  NPV(xk,yk)+  Ak(b-g(xk,yk) ) . 

If  0(Ak)  is  close  enough  to  the  lower  bound,  then 
stop,  return  xk,yk.  Else,  go  to  step  4. 

Step  4.  Let  Ak+1=  max{0,  Ak  + tk (g (xk, yk) -b)  } where 

_ e*(0  (A*)  -0*  ) 

»<7*ll2 


and  0 < e < 2 , 


and  Ak+1 


(All,A^»  * * * > ^tav'  ' ‘ ‘ ' ^jv)  ’ 
Update  Atav  for  all  av,  av=l, 


. .,V  by  using  eg. 4. 4. 2. 


Set  k=k+l  and  go  to  step  1. 

4.4.3  Computational  Results 

The  Subgradient  algorithm  was  programmed  in  FORTRAN  and 
executed  on  IBM  3090/400  with  vector  processing.  The  GBD 
procedure  coded  by  Erenguc,  Tufekci  and  Zappe  was  used  as  a 
subroutine  to  solve  problem  (SA) . For  the  details  of  the 
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subroutine  refer  to  Erenguc  et  al.  [19].  Initially,  the 
lagrangian  multipliers,  X,  are  set  to  zero.  The  heuristic 
solutions  found  at  section  4 . 3 were  used  as  the  lower  bounds 
(LB)  on  0 . A tolerance  of  5%  was  used.  That  is,  the 
algorithm  terminates  when  the  solution  obtained,  upper  bound 
(UB) , is  within  5%  of  the  lower  bound.  The  maximum  number  of 
iterations  was  set  to  20  and  the  17  of  the  50  test  problems 
were  used  to  test  the  efficiency  of  the  procedure.  For  each 
problem  the  average  number  of  GBD  iterations  (i.e.,  the  number 
of  times  the  relaxed  master  problem  was  solved  in  each  call  of 
the  GBD  subroutine)  is  found  to  be  one  and  the  rest  of  the 
computational  results  is  summarized  in  Table  12. 

In  the  table,  column  two  gives  the  heuristic  solutions 
found  for  the  test  problems  which  constitute  the  lower  bounds. 
The  third  column  gives  the  lowest  upper  bound  values  found 
within  computational  limits  for  each  problem  and  column  four 
gives  the  corresponding  lagrangian  multiplier  values.  The 
columns  five,  six,  and  seven  summarize  the  computational 
efficiency  of  the  GBD  subroutine:  Average  number  of  Grinold 
pivots  done  to  solve  a subproblem  is  given  in  column  5,  the 
average  number  of  linear  problems  evaluated  to  solve  the 
relaxed  master  problems  of  GBD  subroutine  is  given  in  column 
6 and  the  average  number  of  simplex  iterations  of  each  linear 
program  is  given  in  column  7.  Note  that  all  of  these  averages 
are  taken  over  the  number  of  subgradient  iterations. 
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The  last  two  columns  give  the  optimality  gap  and  the  CPU 
times  consumed  for  each  problem.  The  following  conclusions 
can  be  drawn  from  the  computational  results: 

The  optimality  gap  increased  with  the  interest  rate.  The 
average  gap  is  9%  with  the  0.01  interest  rate  and  21%  with  the 
0.02  interest  rate.  The  increase  in  the  gap  as  a result  of  an 
increase  in  the  interest  rate  is  due  to  the  fact  that  as  the 
interest  rate  increases,  the  cost  of  completing  an  activity 
earlier  or  later  than  its  optimal  completion  time  costs  more. 

The  computational  results  are  encouraging  in  that  16  out 
of  17  problems  required  less  than  9 seconds  to  obtain  the  best 
upper  bounds.  Only  one  problem  required  12  seconds.  This 
computational  efficiency  solely  depends  on  the  efficiency  of 
the  GBD  procedure. 

4.4.4  Conclusion 

In  this  section,  we  have  presented  a Subgradient 
Optimization  approach  to  RCTCTPSPDC.  Computational  experience 
is  reported  for  17  test  problems.  The  results  demonstrate 
that  with  lower  interest  rates  the  Subgradient  algorithm 
obtains  closer  upper  bounds  to  optimal  solutions  with  modest 
computational  effort. 

4.5.  Concluding  Remarks 

In  this  chapter,  we  have  introduced  and  formulated  the 
resource  constrained  time/cost  tradeoff  project  scheduling 
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problem.  We  have  further  introduced  and  discussed  in  detail 
the  solution  algorithms,  exact  and  heuristic,  for  solving  the 
RCTCTPSPDC.  To  our  knowledge  this  study  constitutes  the  first 
attempt  to  research  the  RCTCTPSPDC.  Computational  experiments 
are  reported  for  50  test  problems  for  the  heuristic  and  17 
test  problems  for  the  subgradient  technique.  The  results 
demonstrate  that  the  heuristic  and  the  subgradient  algorithms 
are  both  efficient  and  practical  methodologies  for  solving 
RCTCTPSPDC. 


Table  12 . Computational  results  with  Subgradient 
Algorithm 


Problem 

** 

LB 

UB 

Lag. 

Multiplier 

Ave. 

Grinold 

solved 

Ave. 

LP 

solved 

Ave. 

Simplex 

Iteration 

% Gap 

CPU 

(sec.) 

i 

30150.17 

33921.14 

5.496 

1 

0 

0 

0.12 

0.6 

5.496 

(14)' 

5.496 

27101.17 

31420.75 

0 

1 

0 

0 

0.15 

0.6 

0 

0 

2 

17213.01 

17393.88 

0 

2 

0 

0 

0.01 

0.01 

0 

(7) 

0 

16705.23 

17218.16 

0 

2 

0 

0 

0.03 

0.01 

0 

0 

3 

21510.37 

23596.43 

0 

1 

0 

0 

0.09 

0.04 

0 

(13) 

0 

18737.10 

22726.99 

0 

1 

0 

0 

0.21 

0.04 

0 

0 

Table  12 — continued 
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Problem 

LB 

UB 

Lag. 

Ave. 

Ave. 

Ave. 

% Gap 

CPU 

Multiplier 

Grinold 

LP 

Simplex 

(sec.) 

solved 

solved 

Iteration 

7 

22398.25 

23008.09 

0 

2 

0 

0 

0.02 

0.10 

0 

(9) 

0 

21849.72 

22739.58 

0 

2 

0 

0 

0.04 

0.10 

0 

0 

8 

7522.07 

8028.10 

0 

4 

0 

0 

0.06 

0.01 

0 

(9) 

0 

6844.85 

7887.39 

0 

3 

0 

0 

0.15 

0.01 

0 

0 

10 

10545.50 

10707.47 

0 

2 

0 

0 

0.01 

0.10 

0 

(10) 

0 

9984.01 

10377.35 

0 

2 

0 

0 

0.03 

0.10 

0 

0 

11 

356.06 

380.15 

0 

2 

0 

0 

0.06 

0.52 

0 

(8) 

0.00083 

2 

0 

0 

0.28 

0.52 

306.53 

371.22 

0.0053 

13 

2220.66 

2665.19 

0 

11 

0 

0 

0.20 

5.55 

0 

(22) 

0 

2001.44 

2726.22 

0 

11 

0 

0 

0.36 

5.55 

0 

0 

17 

2013.98 

2405.75 

0.00019 

8 

1 

10 

0.19 

12.0 

0.0050 

(22) 

0.0022 

1590.44 

2184.39 

0.174 

8 

1 

10 

0.35 

12.0 

0.287 

0.095 

19 

3573.49 

4088.66 

0.265 

7 

0 

0 

0.14 

8.1 

0 

(22) 

0 

2968.62 

3891.35 

0.107 

7 

0 

0 

0.31 

8.2 

0 

0 

Table  12 — continued 
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Problem 

+* 

LB 

UB 

Lag. 

Multiplier 

Ave. 

Grinold 

solved 

Ave. 

LP 

solved 

Ave. 

Simplex 

Iteration 

% Gap 

CPU 

(sec.) 

20 

(22) 

1654.58 

1733.35 

0.013 

0.0010 

0.0020 

7 

0 

0 

0.04 

0.67 

1345.58 

1640.30 

0.0023 

0.0047 

0.0009 

7 

0 

0 

0.21 

0.67 

21 

(22) 

3848.36 

4386.61 

0.000 

0.011 

0.0001 

10 

0 

0 

0.14 

8.0 

3209.34 

4207.73 

0.202 

0.018 

0.020 

10 

0 

0 

0.30 

8.1 

24 

02) 

4612.62 

5559.36 

0.108 

0.031 

0.009 

2 

0 

0 

0.20 

5.9 

3655.90 

5159.14 

0.202 

0.049 

0.016 

2 

0 

0 

0.40 

5.9 

26 

(22) 

4773.47 

5126.06 

0.051 

0.005 

0.002 

4 

0 

0 

0.07 

5.2 

4196.13 

4893.28 

0.047 

0.098 

0.001 

4 

0 

0 

0.16 

5.2 

27 

02) 

4099.27 

4503.51 

0.092 

0.043 

0.001 

5 

0 

0 

0.09 

5.8 

3462.69 

4308.05 

0.076 

0.028 

0.005 

5 

0 

0 

0.24 

5.7 

44 

02) 

4380.30 

4779.19 

0.0031 

0.0511 

0.0166 

4 

0 

0 

0.09 

4.3 

3313.08 

4608.72 

0.000 

0.000 

0.000 

5 

0 

0 

0.38 

4.3 

45 

02) 

4544.57 

4778.84 

0.000 

0.001 

0.004 

10 

0 

0 

0.05 

7.6 

4239.79 

4673.60 

0.029 

0.033 

0.010 

11 

0 

0 

0.10 

7.6 

**  averages  are  per  subgradient  iteration 
* number  of  activities 


CHAPTER  5 

SUMMARY  AND  CONCLUSIONS 


In  this  dissertation,  we  have  suggested  solution 
techniques  for  two  new  project  scheduling  problems:  The 
Resource  Constrained  Project  Scheduling  Problem  with 
Discounted  Cash  flows  and  The  Resource  Constrained  Time/Cost 
tradeoff  Project  Scheduling  Problem  with  Discounted  Cash 
flows.  The  first  problem  involves  scheduling  the  project 
activities  with  cash  inflows  and  outflows,  in  such  a way  that 
the  net  present  value  of  cash  flows  is  maximized  subject  to 
resource  and  precedence  constraints.  A Tabu  search  procedure 
(Tabu-S)  was  proposed  as  a solution  technique  for  this 
problem.  This  procedure  was  modified  to  invoke  a long  term 
memory  function  (Tabu-L) . Both  procedures  were  tested  on  50 
test  problems.  Solutions  obtained  by  these  procedures  were 
compared  to  upper  bounds  obtained  from  a linear  programming 
relaxation  of  RCPSPDC  that  is  strengthened  by  valid  cuts. 
Furthermore,  a comparison  of  these  solutions  to  solutions 
obtained  by  Minimum  Slack  heuristic  was  provided.  In  general, 
it  was  concluded  that  for  small  size  problems  Minimum  Slack 
Heuristic  produced  slightly  better  results  for  problems  with 
activities  that  have  mostly  cash  inflows.  It  was  further 
concluded  that  Tabu  search  is  a viable  approach  for  RCPSPDC  in 
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that  for  all  the  problems  solved,  good  quality  solutions  were 
obtained  with  modest  computational  effort. 

A depth-first  Branch  and  Bound  Algorithm  is  presented  for 
solving  the  RCPSPDC  optimally.  The  algorithm  constitutes  the 
first  effort  to  solve  RCPSPDC  optimally.  Computational 
testing  on  50  problems  demonstrated  that  the  solutions  to 
small  size  problems  were  obtained  within  0.5%  optimality  with 
reasonable  computational  effort.  For  large  size  problems  best 
lower  bounds  to  solutions  were  obtained  within  30  minute  CPU 
time  limit.  In  general,  it  was  concluded  that  the  Branch  and 
Bound  algorithm  can  be  used  effectively  for  solving  moderate 
size  RCPSPDC.  Furthermore,  for  the  test  problems  considered, 
the  results  demonstrated  that  there  is  almost  always  less  than 
1%  optimality  gap  between  the  Tabu-L  solutions  and  the  Branch 
and  Bound  solutions. 

The  second  project  scheduling  problem,  RCTCTPSPDC, 
involves  determining  the  duration  and  scheduling  of  activities 
with  cash  flows  in  such  a way  that  the  resource  and  precedence 
constraints  are  met  and  the  net  present  value  of  the  cash 
flows  is  maximized.  This  problem  was  first  introduced  and 
formulated  in  this  dissertation.  A heuristic  procedure  was 
first  developed  which  schedules  the  activities  according  to 
their  Net  Marginal  Gain  values.  The  computational  testing  on 
the  same  set  of  test  problems  demonstrated  that  it  can  produce 
feasible  solutions  to  not  only  small  size  but  also  to  very 
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complicated  and  larger  project  scheduling  problems  with  small 
computational  effort. 

In  addition,  a Subgradient  Algorithm  which  employs 
Generalized  Benders  Decomposition  technique  as  a subroutine 
was  developed  for  obtaining  an  upper  bound  for  RCTCTPSPDC. 
The  solutions  to  test  problems  were  obtained  which  are  within 
9%  and  21%  optimality  with  0.01  discount  rate  and  0.02 
discount  rate  respectively.  The  procedure  performed 
efficiently  mainly  due  to  the  efficiency  of  the  GBD  algorithm. 

A natural  extension  of  the  work  done  in  this  dissertation 
would  be  to  develop  and  apply  new  methodologies  to  RCPSPDC  and 
RCTCTPSPDC  in  order  to  improve  and  compare  the  efficiency  of 
the  suggested  algorithms  to  be  able  to  solve  large  scale 
problems.  Furthermore,  Tabu  Search  algorithms  can  be 
developed  for  solving  RCTCTPSPDC. 


APPENDIX  A 


THE  FIXED  DEADLINE  ALGORITHM 

Notation: 

<5  : The  upper  bound  on  the  duration  of  the  entire  project. 

<50  : The  briefest  possible  duration  of  the  project. 

V : The  set  of  nodes,  V = { 1,...,N). 

Tj  : The  realization  time  of  event  i. 

r : The  arcs  in  the  project  tree. 

t ' : The  arcs  deleted  from  A to  form  the  tree. 

Di j : Duration  of  activity  (i,j). 

q;  : cash  flow  at  the  realization  of  event  i. 

The  Algorithm: 

Phase  1.  Find  the  early  tree  and  the  early  start  times,  T. 
Set  Tn  = <50.  If  <S0  > <S  no  schedule  can  meet  deadline. 
Otherwise,  calculate  b(  = q{  exp(-aTj)  for  i = 1,2,...,N;  set 
PV  = E bj . 

Phase  2 . 

Step  1.  Finding  an  arc  to  delete  from  the  feasible  tree.  Set 
Wjj  =0  for  (i,j)  e r 1 , then  solve 

If  all  w^  >=0,  the  current  solution  is  optimal.  If  some  wtk 
<0,  go  to  step  2. 
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Step  2.  Finding  a new  arc.  Deleting  arc  (i,k)  from  (V,r) 
forms  subgraphs  (U,L)  and  (W,R)  where  node  1 is  in  U.  If  1 eU 
go  to  step  3;  if  1 eW  go  to  step  4. 

Step  3.  Delay  the  events  in  W.  Compute 

0 = Min  [Tj  - Tj  -Djj  | i eW,  jeU,  (i,j)  er']. 
Suppose  arc  (m,n)  achieves  the  minimum.  Go  to  step  5. 

Step  4.  Advance  the  events  in  W.  Compute 

0 = Max  [TV  +Dn  - Tj.  ! i eU,  j eW,  (i , j ) er]. 
Suppose  arc  (m,n)  achieves  the  maximum.  Go  to  step  5. 

Step  5.  Change  of  schedule.  Let  A = exp(-a0).  For  i eW, 

Tj  <-  Tj  +0 

bj  <-  bj  A 

Delete  (l,k)  from  r and  add  (m,n).  Go  to  step  1. 

The  above  algorithm  finds  the  optimal  solution  to  the  payment 
scheduling  problem  for  any  given  project  deadline. 


APPENDIX  B 


GENERALIZED  BENDERS  DECOMPOSITION 
Generalized  Benders  Decomposition  (GBD)  technique  can  be 
used  for  solving  nonlinear  programming  problems  in  the 
following  form: 

(NP)  Maximize  f(x,y) 
subject  to 

G (x, y)  > 0 
x e X and  y e Y. 

The  vector  y is  called  the  complicated  variables  such  that  if 
y is  temporarily  held  fixed  then  the  problem  (NP)  can  be 
solved  easily  in  X.  The  GBD  procedure  can  be  summarized  as 
follows  [19]: 

Step  1.  For  a given  y e Y,  solve  the  following  optimization 
problem: 

(sub)  Maximizex  f(x,y) 
subject  to 

G(x,y)  > 0 
x e X. 

The  optimal  solution  to  problem  (sub)  provides  a lower  bound 
(LB) on  the  optimal  solution  to  (NP)  problem.  Let  u>0  be  the 
optimal  multiplier  vector  associated  with  the  G(x,y)>0 
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constraints.  It  is  assumed  that  problem  (sub)  is  solvable  by 
a dual  adequate  algorithm.  Let  p=l  and  u1=u,  go  to  step  2. 
Step  2.  Solve  the  following  optimization  problem: 

(RMP)  Maximize  yo 

y e Y>Y0 

subject  to  yo  < maximum  (f(x,y)  + uj  G(x,y) ) j=l,...,p 

x e X 

Let  (y*,y*0)  be  the  optimal  solution  to  (RMP)  . Let  UB=y*o, 
where  UB  is  an  upper  bound  on  the  optimal  solution  to  (NP) 
problem.  Go  to  step  3. 

Step  3.  Set  y=y*  and  solve  problem  (sub)  again.  If  LB>UB-d 
where  d is  a convergence  tolerance,  stop,  the  solution  to 
(sub)  is  the  solution  to  (NP) , otherwise  p=p+l  and  up=u.  Go 
to  step  2 . 

In  his  paper  Geoffrion  has  proved  that  the  above 
procedure  terminates  in  a finite  number  of  iterations  for  any 
given  convergence  tolerance,  when  Y is  a finite  discrete  set 
and  problem  (sub)  has  an  optimal  multiplier  vector. 

The  Applicability  of  GBP  to  TCTPSPDC 
Erenguc,  Tufekci  and  Zappe  in  their  paper  [19]  have 
showed  the  applicability  of  GBD  technique  to  time/cost 
tradeoff  project  scheduling  problem  with  discounted  cash 
flows.  In  what  follows  is  a brief  summary  of  the  paper: 

The  formulation  given  in  section  4.4.2  was  rewritten  in 
a form  that  is  more  amenable  to  solution  with  GBD.  In  this  new 
form,  the  different  patterns  of  durations  of  the  arcs  are 
determined  that  are  leading  in  to  events.  Thus,  the 
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formulation  now  has  new  variables,  y jk,  which  takes  a value  of 
one  if  the  kth-  pattern  of  durations  is  chosen  for  some  event 
and  zero  otherwise.  Next,  variable  transformations  on  the  new 
formulation  were  performed  so  that  when  y is  fixed  the  problem 
becomes  a linear  programming  problem.  Thus,  when  the 
complicated  variables  are  fixed,  the  optimal  set  of 
multipliers  was  generated  by  dual-adequate  algorithm, 
Grinold's  fixed  deadline  algorithm  [33].  A depth-first  branch 
and  bound  algorithm  was  developed  for  solving  the  (RMP) . At 
each  node  of  the  tree  a heuristic  procedure  was  implemented  to 
find  a feasible  solution  to  (RMP) . Computational  testing  on 
project  networks  demonstrated  the  efficiency  of  the 
application. 
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